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Disk resistors between | and 10° milliohms needed for the radio-frequency micropotenti- 


ometer were not available. 
processes for forming thin-film elements. 


Therefore, considerable research was necessary on a wide variety of 
Those methods whereby resistors stable to within 


plus or minus 1 percent were consistently produced are described in detail, and the various 
limitations, fabrication problems, and solutions to these problems are discussed. 


l. Introduction 


With the rapid growth of the electronics field, 
particularly in the past decade, there has been an 
increasing need for accurate voltage standards in the 
microvolt range for radio-receiver sensitivity meas- 
urements, field-intensity determinations, the calibra- 
tion of r-f signal generators, and other applications. 
\t the present time an accuracy of +1 percent has 
become highly desirable. A device developed to 
meet this need, called the r-f micropotentiometer,’:* 
provides an accurate, portable, and rugged low- 
impedance source of r-f voltages from 1 to 10 ° uv for 
use both in the laboratory and in the field. To 
achieve a low-output impedance of 1 to 10° milli- 
ohms, it was necessary to develop low-resistance disk 
elements. Research was conducted on a great many 
different processes for making thin-film resistive 
elements, to determine which consistently produced 
elements having the requisite stability of +1 percent 
over long periods of time and a satisfactorily flat 
frequency characteristic to at least 300 Me. 


The following processes were tried: (a) clamping of | 


resistive disks, (b) low-temperature painting and 
baking, (c) high-temperature firing, (d) evaporation 
and plating, and (e) chemical reduction. Only three 
processes will be described here, that is, high-temper- 
ature firing, evaporation and plating, and clamping. 
Low-temperature paimts proved unsatisfactory, 
whereas the value of the chemical-reduction process 
has not been fully determined. 


2. General Description 


Basically, a micropotentiometer consists of a disk 
resistor element whose r-f impedance is resistive, 
the same as for direct current, and some means to 
indicate the potential across that element. For ex- 
ample, figure 1, a, shows such a resistor terminating 
an essentially lossless coaxial line. The input voltage 
to the line, (1 to 10 v), may be measured accurately 
by means of a calibrated r-f voltmeter. The current 
at the resistive element is 
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Tr Z,, cos Bl+-jZp sin al’ 


where at a given frequency 
[,=current at the resistive element, amperes 
E 
Z 
Z 


.s=r-f voltage input to the line, volts 

i=Iimpedance of the resistive element, ohms 

o=characteristic impedance of the line, ohms 

B=2x/r 

/=electrical length of the line, centimeters 
A=wavelength of the r-f signal in the line, 

centimeters. 


Therefore, 
E’sZ1 
Z, cos Bl+-jZ, sin Bl’ 


Ep=1IprZ, 


where Ee» is volts at the resistive element. At those 
frequencies where / is an odd multiple of 4/4, eq (2) 
reduces simply to 

E'sZ1 


Z. (3) 


ER 


Figures 1, b, and 1, ¢c, illustrate a thermoelement and 
a bolometer, respectively, used to measure the current 
flowing in the resistor, R. For these latter applications 
the current monitor is considered as having no fre- 
quency error even at the highest radio frequency used. 

By using resistances of from 1 milliohm to 1 ohm 
in combination with current levels between 1 and 
100 ma, levels from 1 to 10° uv can be made available. 
The voltage range may be extended to lower or 
higher values. However, the use of higher resist- 
ances may make it necessary to correct for the 
impedance of some external circuits connected across 
the resistance element. 

The resistance of an annular ring is principally a 
function of its physical dimensions and the electrical 
conductivity of the resistive material used. Ex- 
pressed mathematically, 


pn 4) 


2rod 
W. Jackson and L. G. H. Huxley, The solution of transmission-line problems 
by use of the circle diagram of impedance, Trans. IEE 91, No. III (London, 1944) 
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Figure 1. Schematics of three types of radio-frequency 


micro potentiometers 


t, disk resistive element; RF output is obtained across the coaxial connector 
1, Coaxial line; b, thermoelement: and c, thermistor 


where 
R =resistance of the annular ring, ohms 
1/p=mho/cm 
r/T; 
thickness of the resistive film, centimeters 
inner radius of the annular ring, centimeters, 
or inches 
outer radius of the annular ring, centimeters, 
or inches 
p=resistivity of the metal, ohms/cm 
A basic sectional line drawing of a typical resistive 
element and the coaxial electrodes with which it is 
in contact, figure 2, identifies the dimensions r2, r;, 
and d in eq (4). When the resistive film was formed 
by firing conductive paints containing the noble 
metals, the resistance values obtained were some- 
what greater than calculations would indicate because 
of the physical properties of the film so formed. 

At zero frequency the current flowing in a con- 
ductor is uniformly distributed over its entire cross 
section. However, as the frequency is increased, 
there is a greater concentration of the current close 


to the surface, thereby increasing the resistance 


above its d-c value. This phenomenon‘, known as 
“skin effect’, must be largely 
resistive element of the micropotentiometer, where 
the convenience of operating with d-c calibration 
values is desirable 

The depth of r-f current penetration, 6, in a thick, 
solid metal is given by 


l b 
( om™S bm ) : 
where 


5=depth of penetration in meters 
o¢,=conductivity of the metal 
}) =frequency of the applied current in ¢/s 
m= 42 <107‘h/m for nonferrous metals. 


(McGraw-Hill Book 


‘F. E. Terman, Radio engineering handbook, p. 4 
Co., Ine., New York, N. Y., 1943 
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Figure 2 Cross-sectional line drawing of resistive elem 


showing the electrodes and the resistive film 


Considering the thin-film annulus as a section of 
coaxial line with a solid conductor as the propagation 
medium, the transfer impedance, as derived from 
transmission-line theory, is, for all practical purposes, 


given by 

- » 7 a 

36° 9064 —-° 
where 


Z,= transfer impedance, ohms 
R,,=d-c resistance of the annulus, ohms 
d=film thickness, centimeters. 
If at any given frequency d<0.56, the r-f current 
penetration is very nearly complete and the r-f 
impedance differs from the d-c resistance by less 
than 1 percent. 


3. Fabrication Problems and Tests 


There are several problems common to all methods 
described in this report for fabricating low resistances 
from thin films. Because these difficulties impos: 
limitations upon the usefulness of the completed ele- 
ment, and may even render it useless, each of these 
problems will be examined. However, the specifiv 
solutions will be presented later in the discussioa of 
fabrication techniques. 











3.1. Radio-Frequency Characteristics and Stability 


The r-f impedance of a resistor element will be fre- 
quency-sensitive if (a) the coaxial electrodes fail to 
contact the film along uniform concentric rings, that 
is, if the contact areas are irregular and spotty so that 
the width of the annulus between them is not uni- 
form, (b) the thickness of the film in the annulus is 
uneven, or (c) the electrodes overlap some part of th 
resistor while remaining insulated from it. Any 
combivatioa of these conditions may be found in a 
single resistive element. 


found, generally condition (c) also exists. The first 


produces an unsymmetrical distribution of curreot in | 


the anoulus resulting in a series inductive reactance 
whereas the latter, because of the capacitance be- 
tween the film and the electrodes, foreshortens the 
r-f current path thereby reducing the r-f resistance 
Therefore, these effects tend to compensate each 
other. Unsymmetrical distribution of the current 


Where condition (a) 1s | 
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he ring is also caused by uneven film thickness. 
tesistor instability is caused by poor electrical 
act between the film and the electrodes, or re- 
; from the flaking off or migration of film particles 
he annulus. After sufficient skill is developed in 
lying the techniques described, the percentage of 
table elements should be very small. Laboratory 
s during the past 2 vears on resistive elements 
le from clamped carbon disks, fired-on conductive 
nts, and evaporated-and-plated films have estab- 
ed that such elements remain stable to + 1 per- 
t or better 


2. Power Dissipation and Mechanical Ruggedness 
Because the micropotentiometec is primarily in- 
nded for use at levels bet ween l and 10° BV, it is usu- 

ily unnecessary to apply more than 100 ma to any 

esistive element. Experimentally, it was found that 

all types are capable of carrying at least five times 
this current without deleterious effects. Therefore, 
there is no power-dissipation problem. 

Resistive elements made by the evaporation and 
plating technique were more fragile and susceptible 
to mechanical injury than those made by the other 
Mechanical shock or abrasion can break 
the film loose from its supporting base (substrate), 
to which it is bonded only by the keying action of the 
film particles. The stresses applied to the center in 
normal use of all types of elements may cause de- 
terioration of electrical contact to the film, and an 
unstable element, unless the coaxial electrodes and 


processes 


the separating dielectric are firmly joined. 


3.3. Temperature Effect 


When the e-f micropotentiometer is used solely in 
the laboratory where the ambient temperature never 
varies by more than a few degrees, the error in re- 
sistance due to temperature variations will seldom 
exceed 1 percent. However, a micropotentiometer 
used in the field may be subjected to a wide range of 
ambient-temperature conditions. In this case, for 
resistive materials haviag a high-temperature co- 
efficient of resistance, errors of 5 to 6 percent are pos- 
sible. ‘To determine how serious this problem might 
become, 1-ohm disk resistors made from carbon sheet, 
25 ohms per square, evaporated, and plated silver 
film resistors and disks made from Hanovia fired-on 


TABLE 1. 


Mean- 
temperature 
coefficient of 

resistance, 
. 
a: 


Resistive material 


via silver paint No. 38 
porated and plated silver 

ia platinum-gold paint No. ! 

ia platinum No. 14 

m sheet, 25 ohms per square 


the percentage in resistance per deg C 


conductive paints (platinum No. 1, platinum No. 14, 
and silver No. 38) were tested over a temperature 
range from 5° to 50°C. The mean-temperature co- 
efficient of resistance, a,, in this temperature range 
for each of the above materials is presented in table 1 


3.4. Factors Limiting the Resistance Range 


For a given resistive material, a resistance range 
of about 10 to 1 may be obtained merely by using 
different film thicknesses. Solid films of metal less 
than 25 mu thick cannot be relied upon because of a 
tendency to break up.° The metal film may form 
into droplets because of the surface tension, or par- 
ticles may migrate until the surface consists of un- 
connected agglomerates of the metal. Film thick- 
ness may be increased as long as the above current 
penetration requirements are met. 


3.5. Initial Tests and Life Tests 


It has been found, experimentally, that structural 
imperfections affecting the stability and frequency 
characteristics of resistive elements do not ordinarily 
alter the d-c resistance sufficiently to be detected by 
simple d-c measurements. A simple and accurate 
initial test that may be used is to incorporate the 
element in a micropotentiometer and compare the 
voltage drop across its output with an accurately 
known r-f voltage or with the output from other 
tested micropoteatiometers. In figure 3,a the known 
reference voltage is provided by another micropo- 
tentiometer, .\/,. Both it and the unknown, WJ, 
are coanected ia parallel branches to the same r-f 
source. Thus frequency drift either in the receiver 
or the generator has negligible effect. The transfer 
instrument is any r-f detector, usually a receiver, 
having a sensitivity of 1 wv, or better. The r-f cur- 
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FIGURE 3. 
unknown micropotentiometer versus a tested one. 


1, Direct comparison; b, comparison of widely different units using a 
standard attenuator 


Robert G. Picard and O. 8. Duffendack, Studies on the structure of thin 
metallic films by means of the electron microscope, J. Appl. Phys. 14, 201-305 
1943 





rent into each micropotentiometer is adjusted until 
the voltage output is the same for both, as deter- 
mined by the transfer instrument. Then each 
thermoelement is calibrated with direct current to 
determine the r-f currents that flowed in the resistive 
elements, and the output voltages are calculated from 
Ohm’s Law. In matching the output voltages the 
current levels may reach extremes that are difficult 
to measure accurately if the resistances of the micro- 
potentiometers are widely different. For example, 
ambient temperature variations materially reduce 
the accuracy with which a current of 20 percent or 
less of the full rated value flowing in a given thermo- 
element can be measured. A standard waveguide- 
below-cutoff attenuator is placed in series with the 
higher resistance micropotentiometer, .\/,, as shown 
in figure 3,b, to keep the r-f-current levels within 
reasonable limits. The r-f input to each branch is 
adjusted to give full rated current in the thermo- 
element of each micropotentiometer, the attenuator 
in the .V/, branch having been set at some convenient 
low attenuation level. The r-f output of VM, is then 
reduced to match that of 4, by adjusting the stand- 
ard attenuator. This is possible because the heater 
resistance of the thermoelement in the VW, branch 
changes by a negligible amount with change ia the 
heater current. 
The final voltage output of VW, is 


‘uy VoXAA, 
where 
\,=voltage drop across the output of M, for 
full rated current, 


change in attenuation, as a voltage ratio 
1, required to reduce V>) to match 


AA-= 


M,- 


oe 


Differences in the output voltages calculated for | 


each of the micropotentiometers are caused by 
structural defects. 

This same test, repeated at intervals of from 2 to 
4 weeks for a period of 6 months, constitutes a test 
for the quality of fabrication techniques when a proc- 
ess has first been put into operation. From actual 
tests it was found that a better indication of quality 
can be obtained at high radio frequencies. 


4. Fabrication Processes 
4.1. Clamped Elements 


A carbon disk clamped to coaxial electrodes pro- 
vided a very stable type of resistor element. Fig- 
ure 4 is a partially exploded view of such an assembly. 
The disk depicted in figure 5 was cut from a com- 
mercial carbon sheet consisting of a very thin car- 
bon film deposited on a backing of Bakelite. Be- 


cause this card material was available in a number | 


of resistance values, it was important to select one 


that would provide an annular resistive element of | 


reasonable physical dimensions. Fabrication of re- 
sistances smaller than 1 ohm was attempted by 
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Ficure 4. Blown-up view of a clamped-type resistive element 
CPL, brass center electrode; GR, brass electrode; C R, brass 
DD, dielectric; G, brass ground shell; | 

N, brass outer ground shell; R, carbor 


CPN, Brass nut 
clamp ring; CN, brass clamp ring 
metal pins; CP, brass center electrode 
coaxial resistor. 


reducing the width of the annulus, but it was ex- 
tremely difficult to obtain satisfactory results when 
the width was made less than 1 mm. As the 
resistance decreases logarithmically with the de- 
crease in diameter ration, y, this method was not 
the correct approach. The reductioa in y necessary 
to produce an appreciable decrease in resistance 
resulted in such close spacing of the electrodes that 
the shunting capacitance between them became 
rather high. Inasmuch as the resistance per square 
was known, it was more convenient to determine 7 
for a particular resistance, R, from eq (8), the 
equivalent of eq (4), except for rearranged constants 


I 
y=log '( oge78 ) 


where S is the specific resistance of the carbon sheet 
in ohms per square. An inver radius, r;, of % im 
was used to provide enough area for the center mount- 
ing hole, and an adequate contact area for the cen- 
tec electrode. The maximum radius of the disk was 
made \ in. larger than r, to provide sufficient contact 
area for the outer coaxial electrode used. 
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Sketch of a carbon disk for the clam ped-ty pe ele- 
showing the annular resistive area and the painted 
f rings 


A;, Painted ring (outside A», painted ring (inside 


There are, undoubtedly, many ways of cutting 
disks from the carbon-coated card material. The 
method used provided quite satisfactory results, 
whether cutting one or a large number of disks at a 
time. The sheet material was placed face down on 
a smooth, flat surface, and squares about \ in. 
larger than the desired outside diameter of the fin- 
ished disk were individually cut on power shears. 
The squares were stacked, clamped between flat metal 
plates, drilled, and turned to size on a lathe. Air- 


drying silver pamts (Linnick Chemical Co. air-drying 
silver lacquer or DuPont silver conducting paint 


No. 4817) were used for concentric areas to be in 
contact with the electrodes. Either was applied 
with a camel’s-hair brush after masking the annular 
ring acea. High precision in this operation was 
found unnecessary. The drying process was accel- 
erated by the application of heat from an infrared 
lamp for a short time. 

Returniag now to figure 4, parts G, D, CR, CN, 
and CPL are shown assembled to form the core of 
the unit. The contact resistance was to be less and 
remain less than 1 milliohm. It was found impor- 
tant that the surfaces of G and CPL against which 
the disk was clamped be coplanar, and parallel with 
the surface F, the faces of GR, and the surface G of 
part N, and that all were normal to the axis of the 
assembly. Stable and positive contact between the 
electrodes and the disk was obtained with tin- or 
lead-foil washers 2 to 5 mils thick. These washers 
were cut to fit exactly over the silver contact areas 
and were placed betwee. the disk and the electrodes 
GR and CP. In clamping the disk iato place, care 
was taken to prevent turaing of the parts. This was 
accomplished by pinning CP to CPL, and GR. to G. 
(he pins used for this purpose fitted loosely so as 
hot to prevent these parts from resting flat against 
the disk. The shank of CP fitted freely into CPL 
and surface F was undercut slightly at C for the 
same reason. In the figure, the edges of the elec- 


trodes appear quite rounded. All that was required, 
actually, was to break the edges enough to remove 
sharp irregularities that could damage the resistive 
film and prevent the electrodes from seating properly. 

Clamped assemblies were found unsatisfactory for 
low-value resistance units. Nichrome foil, which 
had suitable resistivity, was too hard. It failed to 
form to the electrodes resulting in spotty contact at 
only a few points. Such a unit, because of lack of 
electrical symmetry, exhibited reactance at relatively 
low frequencies. The noble metals that were re- 
quired as thin foils (2.5 yw thick for resistance of 1 
to 10 milliohms) proved to be too soft. Under the 
pressure of contact, they flowed sufficiently to cause a 
deterioration of contact to the point where the con- 
tact resistance became comparable with the foil 
resistance. 


4.2. Fired Paints 


Where stable elements of less than 1 ohm were 
needed, the construction illustrated in figure 6 
produced good results. <A dielectric, D, of glass or 
ceramic, completely covered, except on one end, by a 
very thin coating of fired-on conductive paint was 
soldered to the inner and outer electrodes G and P. 
That portion of the film covering the end of D 
constituted the resistive annular ring. Inasmuch as 
the electrodes were soldered directly to the metal 
film comprising the resistive element, a properly 
soldered unit was free from the contact resistance 
difficulties encountered with the clamped assembly. 
Many conducting paints and pastes of platinum, 
gold, silver, and combinations of these were tried. 
Those listed in table 2 were selected because resistive 
films made from these possessed the requisite stability 
adhered well to their substrates, and had less tend- 
ency to alloy with the solder in the assembling proc- 
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Fioure 6. Cross-sectional view of an element using a glass or 
ceramic sleeve metallized by the fired-paint tech nique. 





TaBLe 2. Approximate data for 


Dilution 
wit! 
turpentine 


Percent 
platinum bright Nome 
platinum 


Hanovia No. 1 
containing gold and 
illoyed 

Hanovia No. 12 
pension 


platinum in sus 5 ! 200 


ter 


using fired paints as resistive elements 


Resistivity Approximate 
From mfg diameter 
data ratios 


Resistance 
obtained 


Firing 
nperature 


Ohmaicm 


0.3 to 0.7 ohms 13 (film 0.1 g thick 


leave door of 18 milliohms 


oven open until 800 


F is reached 


Hanovia No. 4 
puste 
Hanovia 14 
Hanovia No. 38 (sily 


platinum alloy Up to Pyrex or ceramic 200 
do an 
Giold and glass or 000 


ceramic 


Unthinned 
er paste Ww 


The information contained in this table was 
intended to serve merely as a rough guide in the 
design and construction of elements. In practice, 
the resistivity of any of these paints varied as much 
as 100 percent from the values shown because it 
depended upon such factors as the degree of dilution, 
the manner of application, the manner of firing, and 
the condition of the surface of the substrate. There- 
fore, only an approximation of y for a particular 
value of disk resistance could be determined from 
eq (4). 

The first step in assembling a fired element was to 
select the proper substrate for the particular con- 
ductive paint to be used. This was important 
because many paints are intended for firing to a 
particular base material. Should this paint be fired 
on a base for which it was not intended, poor bonding 
to that base could result, and the finished element 
would be mechanically fragile. Because the paints 


ess. 


fired at a ag nce of approximately 1,400° F, 


only a Pyrex glass or a ceramic material capable of 
withstanding this temperature without undue stress- 
ing or distortion could be used. 

Dielectric bushings with end faces perpendicular 
to their axes were cut from glass or ceramic tubing, 
using a thin glass-cutting wheel. One end face was 
polished to provide the smoothest possible surface 
and its edges were slightly rounded. In the case of 
glass, this was easily done by fire-polishing the end 
with the flame of a torch; with ceramic, this was 
much more difficult and the finished surface was 
rather porous. Consequently, commercially avail- 
able glazed ceramic insulators of the style illustrated 
in figure 7, a, were used. With this configuration for 
D, the outer electrode, G, needed modification, as 
shown in figures 7, b, and 8. 

Before applying conductive paint to the insulator, 
the surfaces were thoroughly cleaned with acetone 
to remove all foreign matter. Once cleaned, the insu- 
lator was handled only with tweezers. A single thin 
coating of conducting paint, properly thinned, was 
applied with a small brush to all surfaces except the 
end having the unrounded edges, and care was taken 
to see that it bridged the rounded edges to form a 
continuous film. For applying the paint uniformly, 
the insulator was held vertically in a motor-driven 
jig, figure 9, and the paint applied as it was slowly 
rotated. Difficulty was encountered in coating the 
inside surface of the center hole when the diameter 


was less than %{, in. 


to 





3) to 43 milliohms 0.75 to O85 when 
thinned 15 percent 

10 to 18 milliohms do 

1 milliohm 0.04 
film 


1,300 O.l-mil thick 


All the paints listed in table 2, with one exception, 
were fired in the same manner. After the paint had 
air-dried until it was no longer sticky, the insulator 
was placed in an electric muffle furnace at room 
temperature. The furnace was then slowly heated 
until the firing temperature was reached. When 
firing Hanovia No. 12, the furnace door was left 
open until a temperature of about 800° F was ob- 
tained to allow a copious supply of air for oxidizing 
the solvents in the paint. At this point the door 
was closed and the heating cycle continued until the 
oven reached the maximum temperature of firing, 
1,400° F. When using the other paints listed, it 
was unnecessary to leave the door open at any stage 
Once the furnace arrived at the maximum tempera- 
ture, it was kept at this heat for about 10 minutes 
and then allowed to cool slowly to prevent internal 
stressing of the dielectric during cooling. The 
platinum-coated insulators were then burnished, 
washed with acetone, silver contact rings painted 
on, and refired as before. 

Prior to soldering the metallized insulator into its 
electrodes, the film was burnished with fine steel 
wool, washed with acetone, and a mask applied to 
the annular portion to prevent the solder from flowing 
over it. Sauereisen Cement No. 78 was used. This 
air-drying compound was thinned with water to the 
desired consistency and applied by brush to cover the 
end face and the rounded edges. If these edges had 
been left unmasked, the solder would have flowed 
over them to form wedge-shaped extensions to the 
electrodes, with the result that the r-f resistance at 
high frequencies would have been more than its d< 
value. After the masking had dried, the exposed 
areas of the film were burnished and washed again 
and Diveo flux No. 300 immediately applied. Any 
delay in applying this flux may result in contamina- 
tion of the surface. 

The next step was to prepare the electrodes for the 
assembling operation. They were placed into a 
duralumin jig, as shown in figures 10, a and b, 
and 11, to position them so that the insulator could 
be dropped into place and soldered in. In addition 
to its function of holding the pieces, the jig conveyed 
heat by conduction to both members simultaneously 
The pieces were heated only enough to cause the 
Diveo 233 solder to melt and flow smoothly. Then 
the surfaces indicated in figures 6 and 7, b, were 
liberally tinned with this solder and the metallized 
insulator dropped into place between the electrodes 
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parts, including the brass outer electrode, the metallized 
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UkE 12. Cross-sectional view of a resistive element using 


a Kovar-glass seald terminal for the electrodes and supporting 
substrate for a resistive film 
Face A is metalized to form the resistive annulus 
and gradually worked about until it was properly 
seated. The fit between the parts was a loose one 
to permit the flow of solder between them to form 
a solid electrical and mechanical bond over the whole 
circumference of each junction. Very low-temper- 
ature melting solder was not found satisfactory as 
it produced a poor mechanical bond. While enough 
solder was used to insure good bonding completely 
around the circumferences of the electrodes, the 
excessive use of solder often resulted in shorting of 
the element, the excess flowing through the unit com- 
pletely and bridging the dielectric on the under side. 

When thinning a conducting paint with turpentine, 
only the amount to be used in a given day was pre- 
pared. Where large quantities were prepared at one 
time and left in storage, the solvents gradually evap- 
orated resulting in a more concentrated mixture 
which had higher conductivity. Moreover, as the 
concentration increased, a single coating produced 
a thicker film. Consequently, elements made at 
different times from the same mixture had widely 
different resistance values. 


4.3. Evaporated and Plated Elements 


Resistive elements of approximately 1 milliohm 
were also constructed by another method, in which 
the resistive film was deposited by the evaporation ° 
and plating of pure silver on a flat surface consisting 
of the coaxial electrodes and the separating dielectric. 


In one type, using this technique, Kovar-glass 
sealed terminals were used, as illustrated in figure 12. 
The proper ratio of diameters, y, for the annular 
ring formed by the glass surface between the elec- 
trodes was determined from eq (4); but, as the ratio 
varied within rather narrow limits for the different 
type terminals available, the over-all physical size 
of the terminal became the prime consideration. 
Only those having a face reasonably flat and free 
from sharp irregularities or breaks, especially at the 
lines of juncture between glass and metal, were 
used. The terminals selected had all metal surfaces 
tinned and had an approximate y of 2:1 and an 
over-all diameter of about \ in. with a center elec- 
trode, which was tubular. Before soldering one 
into its brass shell, the center electrode was shaped 
at both ends into fingered female pins similar to a 
type N. 

*H. W. Weinhart and H. G. Wehe, Metal electrodes deposited on quartz 


crystals by the evaporation process; R. A. Heising, Quartz crystals for electrical 
circuits, p. 333-355 (D. Van Nostrand, New York, N. Y., 1946 


Better elements resulted when the evaporated and 
plated silver film was deposited on_ silver-plated 
electrodes. Each step necessary in depositing the 
evaporated and plated film is listed in the following 
text. To silver plate the electrodes: 

a) The tinning was dissolved from 

with HCI. 

b) The exposed metal surfaces were scrubbed 
under tap water using a fine abrasive, and 
rinsed. 
copper flash was applied immediately in a 
copper-cyanide plating bath (a necessary 
step since unprotected Kovar oxidized very 
quickly if left exposed to the air). 

After rinsing in tap water, a silver strike was 
applied in a standard silver-cyanide strike 
solution and this was followed by plating 

2 to 5 ten-thousandths of an inch of silver 
in a standard cyanide plating bath. 

After rinsing once more in tap water, the 
surfaces were cleaned with acetone and dried 
Then the resistive film was formed as follows 

A thin silver film approximately 0.1 y» thick 
was evaporated upon the entire surface. 

The silver was built up to a total thickness of 
no more than 2.5 w by plating in a standard 
high-concentration plating bath. 

Steps e, f, and g were performed as quickly as pos- 

sible. A delay between these successive steps of as 

little as 1 hour was enough for contamination of the 

surfaces to occur, resulting in poor electrical contact 

| between the film and the electrodes. A final protec- 
tive coating of lacquer, Krylon, was applied to the 
surface to prevent mechanical injury. 

The second type of unit was similar to that 
described in the section devoted to the fired-paint 
technique. After the metallized insulator was 
soldered between the electrodes, the conductive film 
on the end of the button was ground away, and the 
électrodes and the insulator were ground until the 
face that they presented was flat and free from flaws, 

| particularly where the metal and dielectric joined 
Once again, all exposed metal surfaces were silver 
plated before evaporating upon them. However, 
when silver was plated directly upon the electrodes, 
it did not cover the thin ring of solder on the surface 
at the lines of juncture between the dielectric and 
metal. All processing steps required in depositing 
the evaporated and plated film on this type unit are 
listed below. To cover the electrodes with a silver 
plate: 

h) The metal surfaces were scrubbed under tap 

water using a fine abrasive and then rinsed 

i) The unit was dipped into 20-percent HC! and 

rinsed again in tap water. 

j) This was followed by a dip into a hot caustic 

solution (2 percent of NaOH per gallon of 
water at 150° F). 

k) The metal surfaces were then scrubbed and 

rinsed in tap water. 

Upon completion of step k, the remaining steps, 
namely, copper flashing, silver flashing and silver- 
plating of the electrodes, the evaporation of the silver 


the Kovar 
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m and the final plating were the same as steps ¢ 
rough g. Aluminum, the only other metal tried in 
is process, gave unsatisfactory results because of | 
wor adherence to the substrata. 


Conclusion and Implications for Further 


Work 


The processes described above are by no means 

ie only ones by which resistive elements may be 

bricated. However, there is. sufficient experi- 
ental data on resistive elements made by these 

iethods to prove their dependability. A chemical 
eduction process, commonly known as the silver- 
muror process, is presently under consideration. 
The steps in the procedure are much the same as 
hose outlined in the last section, with the exception 
that the chemical reduction of silver replaces the 
final steps of evaporation and plating. 

From preliminary tests on a recently completed 
group of 1- to 5-ohm elements made with fired-on 
paints, it appears that this may be a simpler method 
than the clamped carbon type for fabricating high- 


resistance elements (1 to 10 ohms). Further tests 
are necessary to determine their stability. 

No attempt was made to reproduce resistance 
values within any specific tolerance. The carbon 
sheet material used in the clamped assemblies has a 
specified tolerance of +10 percent, and consequently 
no element of this construction can be expected to 
be reproduced closer than this. The process that 
sockahiy would lend itself best to close resistance 
tolerances is the one using the fired-on paints. It is 
conceivable that by controlling the degree of dilution 
and by using some other means of applying the 
paint, such as spraying or dipping, resistance values 
could be reproduced to within 5 to 10 percent. 


Credit should be given to Vernon A. Lamb and 
William H. Metzger, Jr., for their valuable coopera- 
tion and technical assistance in investigating the 
various chemical problems encountered in making 
thin metallic film. 


WasHINnGTon, May 1953, 








Properties of Some Masonry Cement 
D. N. Evans, A. Litvin, A. C. Figlia, and R. L. Blaine 


Seventeen brands of masonry cement were tested for soundness, strength, consistency, 
fineness, time of setting, shrinkage, water repellency, water retention, autoclave expansion, 


air entrainment, and resistance to freezing and thawing 
Autoclave tests indicated potential unsoundness in some of the 
Air entrained in the mortar was effective in increasing the ability to withstand 
The shrinkage of these 


by weight and by volume 
cements 
freezing and thawing 


Mortars were proportioned both 


cements differed greatly The results of 


tests indicated also that the properties of these masonry cements were appreciably different 


from those previously reported in 1934 


1. Introduction 


\lasonry cements are used in building construction 
to a greater extent today than ever before. Possible 
reasons for the popularity are their general avail- 
ability, the convenience of proportioning and mixing 
masonry mortars, and their desirable physical prop- 
erties. Evidence of appreciable changes in the 
properties of these cements since the 1934 National 
Bureau of Standards survey of these materials [1]* 
has been presented by Watstein and Seese [2] and 
also has been evident in the acceptance-testing-work 
of the Bureau laboratories. In addition to the de- 
velopment of the materials, new methods and ideas 
of testing have been devised since the last major 
revision of the Federal Specification for Masonry 
Cements. It seemed desirable, therefore, to apply 
the knowledge of these new test methods and ideas 
of testing in making a survey of the properties 
of modern masonry cements before revising the 
specifications. 


2. Scope of Investigation 


Tests were made on 17 brands of masonry cement, 
purchased for the most part from building supply 
dealers east of the Mississippi River. 

The following tests were made as required in 
Federal Specification SS-C-181b for masonry ce- 
ment: Normal consistency, time of setting, soundness, 
strength, water retention, and water repellency. 

In addition to these specified tests, determinations 
were also made of the fineness by the No. 325-sieve 
and air-permeability methods, soundness of neat 
cement pats in steam and of neat cement and mortar 
bars in the autoclave test, air entrainment, dura- 
ability in freezing and thawing tests, and volume 
change of mortar bars and reinforced mortar bars. 
The freezing and thawing tests were made on mortar 
bars proportioned with 1 part of cement to 3 parts 
of blended sand, by volume. The other mortar tests 
were made on mortars proportioned both by volume 
and by weight. 


gures in brackets indicate the literature references at the end of this paper. 
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3. Materials Tested 


Cements. The 17 brands of masonry cement 
tested were purchased in 1947 at a time when cement 
was scarce, so these cements probably had no ware- 
house deterioration, and represented the product as 
manufactured. An attempt was made to get as 
many different types of cements as possible. The 
present selection of brands is not, however, repre- 
sentative of the quantities of the various types, as 
the authors have been informed that approximately 
90 percent of the masonry cement sold at present ts 
of the portland-cement-base type. Of the cements 
tested, Nos. 6 and 12 are believed to be largely 
natural cement; No. 7 is believed to contain some 
natural cement; No. 8 is believed to be a slag cement; 
and No. 15, a mixture of portland cement and lime. 
The others are believed to be mixtures of portland 
cement and limestone or raw mix. 

Sand. The sands used were standard Ottawa silica 
sands, both “20-30” and “graded”. The 20-30 
Ottawa sand was used in the air-entrainment mortars 
(1 to 4) and in the mortars for compressive-strength 
tests (1 to 3, by weight). Mortars proportioned by 
volume were prepared with | part of cement and 3 
parts of blended sand (one-half 20-30 Ottawa sand 
and one-half graded Ottawa sand). 


4. Method of Testing 


(a) Fineness on a No. 200 sieve, time of setting, 
soundness of pats in air and water, compressive 
strength of cubes (1 to 3, by weight) water retention, 
and water repellency were determined by the methods 
given in Federal Specifications SS-C—181b. 

(b) The following tests were made according to 
the methods given in Federal Specification SS—C- 158, 
with the exceptions as noted: Fineness on No. 325 
sieve and by air-permeability method, soundness of 
pats in steam and of neat cement bars and mortar 
bars in the autoclave,' and air entrainment.’ 

(c) The following tests were made of the mortars 
proportioned by volume: Autoclave expansion at 7 


Specimens were removed from the molds at 48 hr and soundness determina- 

tions were made at 48 hr rather than 24 hr 

? The mortar in the air-entrainment test 
instead of a rod. 


was tamped with al ty \¢ in. tamper 





days, compressive strength, freezing and thawing. 
water retention, water repellency, and volume 
change of mortar bars and reinforced mortar bars. 

The weights of the materials for the 1 to 3 by 
volume mortars were calculated by assuming the 
weight of a bag of cement to be the weight per cubic 
foot and assuming the dry weight of sand in 3 ft 
of loose damp sand as being 240 lb. This proportion 
was selected because it is the usual recommendation 
in mortar specifications to mix 1 sack of masonry 
cement with ft* of sand for masonry mortar. 
The sand used was a blend of one-half 20-30 and 
one-half graded Ottawa sand, with a fineness modulus 
of approximately 2.40. This sand is also the same 
as that required in the standard method of test for 
measuring mortar making properties of fine aggregate 
ASTM designation C87—47. This blend more nearly 
approximates masonry sand than does 20-30 Ottawa 
sand specified for the compressive strength test in 
Federal Specification SS—C-—18 1b. 

The mortars were mixed in a model C-10 Hobart 
mixer with a special paddle constructed with side 
strips of bronze attached to a center strip. The side 


9 
” 


strips were covered with sections of rubber tubing of 


such length that the ends brushed the side of the 
3-qt bowl in which the mortar was mixed. The 
cement was added to the water in the mixing bowl 
and mixed for 30 sec. About half the sand was added 
and mixed for 30 sec. The rest of the sand was 
then added and mixing was continued for 60 sec. 
The mortar was then allowed to stand for 90 sec, 
with the bowl and mortar covered with a damp cloth, 
and then was remixed for an additional 60 sec, and 
the cubes were made in the regular manner. 

It should be noted that the mortars used in this 
investigation are not the same as the mortars used in 
the investigation by Rogers and Blaine [1], and hence 
the results are not directly comparable. Of the 3 


mortars used in the earlier investigation, 1 mortar | 


was proportioned by weight, 1 part of cement to 3 
parts of 20-30 Ottawa sand, the othe 
portioned 1 part of cement to 3 parts of mixed sand, 
1 by weight and the other by rodded volume. The 
mixed sand was composed of 2 parts of 20-30 Ottawa 
sand to 1 part of pit-run Ottawa sand, the mixed 
sand having a fineness modulus of 2.62. The sand 


and cement were hand mixed in a metal vessel, water | 


was added and the mixing continued for 5 min. 
After standing for 14% min, they were again stirred 
for % min. 

In the tests described herein, cubes were stored in 
the moist cabinet for 7 days, after which specimens 
for test at 28 days and 1 year were immersed in water 
at 70° F. Specimens to be tested after 1 year in 
air were stored in the laboratory air at 70° F and 
50-percent relative humidity after the first 7 days. 
Six cubes were tested at each age after each condition 
of storage. 

Specimens for the freezing and thawing test were 
mortar prisms 2 by 2 by 12 in. After 48 hr they 
were removed from the molds and stored in water 
for 12 days, then in laboratory air for 14 days. At 
the age of 28 days they were placed in water in 24 by 


2 were pro- | 


| 
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24% by 13-in. noncorroding metal containers and 
stored at 40° F for 48 hr. Then the containers hol:- 
ing the specimens were placed in rapidly circulating 
ethylene glycol at —20° F for 2 hr, then in the 40° F 
bath for 2 hr. Two freezing cycles were made each 
day. Specimens were stored at 40° F overnight and 
over weekends. After each five freezing cycles, each 
bar was weighed and sonic measurements [3, 4] were 
made. Measurements were continued to 100 cycles, 
except for bars broken prior to 100 cycles. Only 16 
cements were used as there was not enough of cement 
2 for this test. 

(d) Two sets of mortar prisms 1 by 1 by 11 in 
were made of the mortar of each cement. Reference 
points were placed in the ends of each of one set of 
prisms for length measurement, leaving an effective 
length of 10 in. between the inner ends of the points 
Cold-rolled steel rods 4 by 11% in., threaded the full 
length except for rounded ends, were imbedded in the 
other set of prisms, the rounded ends serving as ref- 
erence points for length change. Molds were re- 
moved at 48 hr, and after an additional 5 days of 
storage in the moist closet, the prisms were placed 
in water at 70° F for 21 days, and then in air at 70° 
F and 50-percent relative humidity. Length meas- 
urements were made at 7 and 28 days and at frequent 
intervals thereafter until apparent equilibrium was 
reached. Examination of the reinforced prisms was 
made for evidence of cracking. 


5. Results of Tests 


Data on weight per bag, density, fineness in terms 
of percentage of weight of cement retained on No 
200 and No. 325 sieves, specific surface as determined 
by air permeability, normal consistency, initial time 
of set, and autoclave expansion of both neat cement 
and mortar are presented in table 1 for all 17 cements. 


TABLE 1. Weight, density, fineness, normal consistency, setting 
time, and linear autoclave expansion of 17 masonry cements 
tested 


Linear autoclave 
expansion 


Weight re- 
tained on 
sieve No 


Nor- 
mal 

con- 
sis- 


tency * 


Fine- 
ness 
(air 
perme- 
ability 


Spe- 

Ce- Bag cific 

ment weight grav- 
ity 


Initial 
~~ 


ting 
time Neat 


NX 4 
cement | Mortar 


325 


b 
70 
70 
75 
7 
70 
65 
70 
sO) 
70 
75 


Ar/min 
4:15 
4:15 
3:20 
3:40 


0.10 
07 
14 
O4 
. 05 


0.11 
.16 
oO 
| 


PrPprT 


2 
4 
03 

. 


0 
. 18 
.08 
7 
.02 


PPP 


.6 
.10 
.0% 
.07 
.o 


70 
62. § 
70 
70 


62 


PRP 


ml 03 
. 


70 


§8 


470 


ne 





* Ratio of weight of water to weight of dry cement times 100 required for pene- 
tration of 10 mm in the Vicat test. 
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BLE 2 Air entrainment, water retention, water re pe llency, 
t-hour absorption and durability and weight loss in freezing 
1 thawing 


Water 
ibsorption 
at 24 hr 


Freezing and 


thawing 


l 1:3 Dura- 
by by bility 
vol - . vol . vol factor 
ume ; "wm [4 


, Weight 


loss 100 
evcles 


1 at 30 cycles, disintegrated at 7 
ht loss at 65 cycles 
lat D eveles 
iat 5 cycles 
rated at 10 cycles 


Compressive at engths and water require ments 


mortars made with Ottawa silica sand 


Shrink 
age of 
Water neat 
re cement 
quire- bars 


vear 
~)| ments 


water 


» 1.73 2,970 49 130 
46 > l 1, 95 2, 010 44 . 116 


f weight of water to weight of dry cement times 100. (Water required 


il flow of mortars). 


Values for air entrainment, water retention, water 
repellency, and water absorption determined on 
mortars proportioned by weight and volume, as well 
as durability factors and weight loss resulting from 
freezing and thawing are given in table 2 for the same 
cements. 

Additional data comprising compressive strength 
at several ages for mortars proportioned 1:3 by 
weight and by volume, together with their water 





Shrinkage , Percent 





~ 
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Time , Days 


Ficure 1. Shrinkage of motar bars (3:1 by volume). 


Stored in air at 70° F. and %-percent relative humidity, Initial length mea 


sured at age of 21 days upon removal from wet storage 


requirements, are given in table 3. This table also 
presents data on shrinkage of neat-cement prisms, 
as well as the water requirement. A selection of the 
data on shrinkage is presented graphically in figure 1. 


6. Discussion 


We ight of Cements, 


The weight per bag ranged 
from 62.5 to 80 lb, and for 11 of the 17 masonry 


cements was 70 lb. In the investigation by Rogers 
and Blaine [1], the weights of 41 masonry cements, 
as rodded into a box, ranged from 39.7 to 89.9 lb/ft 

the mean weight being 64.7 lb. No direct compar- 
ison can be made, as the rodded weight per cubic 
foot was not determined in this study, and the specific 
gravity was not determined in the 1934 investigation. 

Specific gravity. Specific gravities of the cements 
ranged from 2.72 to 2.94, apprecie bly ‘ower than the 
specific gravity of portland cement, which is approx- 
imately 3.15. 

Sieve fineness. Residues on the No. 200 sieve 
ranged from 0.0 to 25.0 percent by weight. Only 
two cements failed to meet the requirement of Fed- 
eral Specification SS-C-181b that not more than 12 
percent of the cement shall be retained on the No. 
200 sieve. Three cements had residues of 15 percent 
or more on the No. 325 sieve. In the 1934 investiga- 
tion the residue on the No. 325 sieve was reported as: 
7 having less then 10 percent, 13 having 10 to 15 per- 
cent, 14 having 15 to 25 percent, 2 having more than 
40 percent, and 3 not reported. Determination of 
the air-permeability fineness was not made in 1934, 





but from a comparison of the sieve residues, it may 
be assumed that the present cements are finer. 

Normal consistency. The water requirement for 
normal consistency ranged from 23.6 to 37 percent, 
and only 3 cements required more than 30 percent. 
The water requirement for normal consistency is the 
ratio of the weight of water to the weight of cement 
expressed as percentage, for the water-cement paste 
to produce a penetratioa of 10 mm by the Vicat 
needle 

Air entrainment. The amount of entrained air as 
determined by the method for portland cement, 
using | part of cement to 4 parts of 20-30 Ottawa 
sand by weight, ranged from 3.7 to 29.0 percent. 
four of the cements had air-entrainment values of 
less than 12 percent, 9 from 12 to 21 percent, and 4, 
more than 21 percent 

Water retention. Of the mortars proportioned by 
weight, only | had a water retention of less than 65 
percent, and of the mortars proportioned by volume, 
t had water retentions of less than 65 percent. A 
majority of the mortars proportioned by weight 
showed slightly higher water retention than the 
mortars proportioned by volume, although there were 
exceptions. 

Aur-permeability fineness. The fineness of the 
cements as shown by the air-permeability method 
ranged from 5,360 to 9,730 em*/g, with an average of 
about 6,780. Masonry cements as a group are much 
finer than portland cement (about 2,700 to 3,500 
cm’/g by the same method). The air-permeability 
fineness of the masonry cements does not appear to 
be related to the sieve residue or to other properties, 
such as strength, water requirements, water retention, 
or water repellancy. 

Time of setting 
3 br 20 min to 7 hr 30 min 
final set in less than 24 hr. 

Soundness. The test for soundness of pats in air, 
water, and steam did not indicate any evidence of 
unsoundness of any of the cements. However, it 
should be pointed out that these soundness tests for 
neat-cement pats are not sufficiently severe to dis- 
tinguish between any of the cements tested. On the 
other hand, the autoclave expansion of neat-cement 
prisms, tested at 48 hr, ranged from 0.03 to 6.2 per- 
cent, with 3 cements having expansion greater than 
0.50 percent, 7 having expansion less than 0.1 percent, 
and 7 having expansion between 0.1 and 0.2 percent. 
The autoclave expansions of the mortar bars propor- 
tioned by vloume and tested at the age of 7 days, 
though lower, fell for the most part in the same order 
as those of the corresponding neat cements. 

In the investigation of white-coat plaster failure 
by Wells, Clarke, Newman, and Bishop [5], it was 
shown that the failures were caused by the delayed 
hydration and resulting expansion of the magnesia 
in the dolomitiec lime. It was also shown by Wells, 
Clarke, and Levin [6], in their study of the expansive 
characteristics of hydrated limes and the develop- 
ment of a test for soundness, that expansion of the 
autoclaved bars was dependent not only on the 
amount of unhydrated magnesia but also on the 


The time of initial set varied from 
All cements reached 





strength of the test specimens. The limits suggest 

by these authors for the autoclave expansion of 1 

cement to lime (by weight) prisms was 1 percent 
In the present investigation of masonry cemet 

only 1 of the 17 cements had an autoclave expansi: 
of more than 1 percent. 

The causes of lack of bond and appearance 
spaces between the brick and mortar bed have bx 
discussed at considerable length by Palmer 
Davis and Troxell [8], and others. A recent paper 
by J. W. McBurney, “Cracking in masonry caused 
by expansion of mortar” [9], questions the assump- 
tion that cracking is necessarily caused by shrinkag: 
He shows that the openings between bricks and 
mortar usually ascribed to mortar shrinkage result 
actually from the mason’s techniques in “cutting 
joints. He reports damaging expansions resulting 
from delayed hydration of free magnesium oxide in 
masonry mortar. 


Water repellency. The ratio of 1-hr absorption 


to 24-hr absorption was less than 0.9 in all cases, and 
ranged from 0.33 to 0.64 for mortars proportioned by 


weight, and from 0.32 
ry 
I 


to 0.76 for mortars propor- 
tioned by he water absorbed at 24 hr 
ranged from 2.5 to 11.1 percent for mortars pro- 
portioned by weight, and from 5.3 to 10.7 percent 
for mortars proportioned by volume. The rate of 
absorption and total water absorbed usually was 
lower for mortars proportioned by weight, although 
there are some exceptions. The results of the water- 
repellency test in this investigation are in marked 
contrast to those reported by Rogers and Blaine 
They tested cements for absorption in mortars X 
(1 cement to 3 parts 20-30 silica sand by weight) and 
mortar Y (1 part cement to 1.5 parts 20-30 silica 
sand and 1.5 parts pit-run silica sand). Of th 
mortars proportioned by weight, one crumbled when 
placed in water and 20 reached 90 percent of the 72-hr 
absorption in | hr or less. Of the mortars propor- 
tioned by volume, 2 crumbled when placed in water 
and 26 reached 90 percent of the 72-hr absorption in 
1 hr or This high rate of absorption was 
related by them to the failure in 10 cycles or less in 
the freezing and thawing test and resulted in the 
development of the water-repellency test that was 
included in Federal Specification SS—-C-18la. They 
also noted that cements having water-repellent addi- 
tions produced mortar having entrained air, and 
that most of these were resistant to freezing and 
thawing. It was also noted that with one exception 
all specimens having low strengths (less than 500 
psi at 28 days) showed poor resistance to freezing 
and thawing, the one exception being a mortar 
having a high percentage of entrained air and a low 
rate of absorption. 

Freezing and thawing. Determinations of the loss 
of weight and of the sonic modulus were made after 
each 5 cycles of freezing and thawing. Mortar 
prisms of cement 15 showed a failure at 5 cycles 
(decrease of 30 percent of the modulus of elasticity 
as shown by sonic measurement) [4] and were s0 
badly crumbled at 10 cycles that no further tests 
could be made. Mortar prisms of cement 8 had a 


volume. 


less. 
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ight loss at 100 cycles of only 9.4 percent. How- 
er, the durability factor was only 9, and the sonic 
asurement indicated failure at 20 cycles. Mortar 
sms of cement 3 had a weight loss of 17.2 percent 
65 eveles and were broken at 70 cycles, and failure 
$ indicated at 30 cycles The three cements that 
led in the freezing and thawing test all had air- 
tertainment values of less than 8 percent. The 
sts were discontinued at 100 cycles. None of the 
ments used in this study had the high rate of ab- 
wption (90 percent of 24-hr absorption in 1 hr) or 
low strength (less than 500 psi at 28 days) that 
re noted in several of the cements tested in 1934. 
(he freezing and thawing test in this investigation 
hows a relation between entrained air and resistance 
» freezing and thawing but does not appear to show 
. good relation between rate of absorption or weight 
iss and resistance to freezing and thawing. Possibly 
. better relation of the tests would have been shown 
had the freezing and thawing tests been carried 
further. The freezing and thawing test is a rather 
severe shock test and may not have the same effect on 
mortar prisms cast in nonabsorbent molds that natural 
eathering has on mortar placed between bricks. 
(om pressive stre ngth. Of the 17 cements tested, 
only cement 6 would be classed as a type I masonry 
cement (not less than 250 and 500 psi at 7 and 28 
days, respectively When the mortar cubes were 
proportioned by volume both cements 6 and 15 failed 
to meet the requirements (500 and 1,000 psi at 7 and 
28 days, respectively) for type Il masonry cement. 
For mortars proportioned by weight, the 7-day 
strength ranged from 370 to 1,980 psi, the 28-day 
strength from 800 to 2,810 psi; and the increase be- 
tween 7 and 28 days ranged from 190 to 1,370 psi, 
or from 13 to 116 percent. The strength of the mor- 
tars proportioned by volume ranged from 310 to 
1,830 psi at 7 days and from 690 to 3,160 psi at 28 
days, with the increase in strength between the two 
ages ranging from 170 to 1,470 psi, or 14 to 122 per- 
cent. For the 12 cements that are believed to be 
portland cement and limestone, the mortars propor- 
tioned by weight had strengths ranging from 1,000 
to 1,830 psi at 7 days and from 1,280 to 2,870 psi at 
28 days, the increase in strength being from 190 to 
| ,040 psi, or from 13 to 90 percent. The compressive 
strength of specimens of 7 of the cements stored for 
vear in air appears to be appreciably higher than 
the corresponding specimens stored in water for 1 
vear. It is known that the degree of saturation of 
a specimen at the time of breaking affects the 
strength. Although the air-cured specimens were 
soaked before breaking it is not known whether the 
degree of saturation was the same in both sets. 
There is also the possibility that some specimens 
that were high in lime or limestone may have been 
weakened by leaching due to prolonzgued water stor- 
age. In the investigation in 1934, no attempt was 
made to bring the strength specimens stored in air 
and water to the same degree of saturation before 
breaking. The conclusion in the 1934 investigation 
vas that air storage usually gave higher compressive 
strengths than water storage. 


Water requirement. The water (expressed as per- 
centage of dry weight of cement) required for a flow 
of 100 to 115 is shown in table 3. The mortars pro- 
portioned by volume appear to require slightly less 
water than the corresponding mortars proportioned 
by weight. The difference is not large enough to 
have much significance due to the wide range of flows 
used (from 100 to 115 The percentage of water 
required for normal flow of the masonry-cement mor- 
tars is about the same as that required for portland 
cement. 

Shrinkage. The unreinforced mortar prisms had 
shrinkages ranging from 0.073 to 0.155 percent in 60 
days. About three-fourths of this shrinkage occurred 
in the first week, with little change after 28 days 
The change of length (shrinkage) against time is 
presented graphically in figure 1 for mortars made 
from cements, which represent the extremes of those 
tested. Note that cements 11 and 15 are alike in 
having high shrinkage but differ in rate of shrinkage 
at both early and late ages. Attention should be 
called to the probability that “The volumetric 
changes of a mortar cast in a non-absorbent mold and 
stored in air is no criterion to the volumetric changes 
that may be expected to occur in the same kind of 
mortar when used as the joint material in brickwork,”’ 
a conclusion reported by Davis and Troxell (8) 

The reinforced mortar prisms showed a shrinkage 
of 0.02 to 0.05 percent as evidenced by compression 
of the reinforcing rods. With the threaded rods, the 
cracks formed were for the most part extremely small 
and difficult to see. This is in contrast to the work 
reported by Watstein and Seese (2), who used regular 
reinforcing rods. Presumably, the regular rods with 
deformations at greater intervals caused larger and 
more widely spaced cracks. The study by Watstein 
and Seese also indicated that for masonry cements of 
comparable compositions, the shrinkage increased 
with the amount of water used. Their work also 
indicated that the resistance to cracking was in- 
creased with increased strength of the mortars and 
was lowered with increased shrinkage of the mortars. 
Their study included the bond between bricks and 
mortars and indicated that both high compressive 
strength and good water retention are important in 
producing high binding strength in brick assemblages 
Adequate compressive strength are required in ma- 
sonry mortars to insure load-bearing capacity, and a 
balance between strength and water retention is 
necessary to obtain a good bond between the brick 
and mortar. 

General. 
or more of a number of different proportions. 


Masonry cements may be made of one 
For a 
cement of a given composition, the relation between 
different properties, such as fineness and strength, 
water requirement, and water retention, can be 


shown. For a large group of masonry cements of 
different materials and various amounts, the rela- 
tionship becomes obscured by the properties of the 
different materials. 

This study indicates that mortars proportioned by 
weight did not always have higher strength, water 
retention, and water repellency than the correspond- 





ing mortars proportioned by volume. Both mortars 
made with cement 8 had the same weight of cement, 
and the higher strength of the mortars proportioned 
by volume is largely due to the difference in sand. 
With the other factors the same, the blended sand 
tends to produce stronger mortars than the 20-30 
sand, as the smaller voids permit a more dense struc- 
ture (1). Some of the differences between the mor- 
tars proportioned by weight and by volume are not 
of sufficient magnitude to be considered significant. 

Some of the differences that may be noted between 
the masonry cements tested in 1947 and those tested 
in 1934 are: a marked tendency toward higher 
strength and slower rates of absorption, waterproof- 
ing agents have probably been largely replaced by 
air-entraining agents, it is probable that ihe majority 
of masonry cements now being made contain air- 
entraining agents, and most of the cements tested 
showed a good resistance to freezing and thawing. 


7. Summary 


Tests were made of 17 brands of masonry cements 
representing materials of different types. The 
results of tests were as follows: 

1. The masonry cements tested had weight-per-bag 
values ranging from 62% to 80 Ib. 

2. The specific gravity of the cements ranged from 
2.72 to 2.94 percent. 

3. The residue retained on a No. 200 sieve ranged 
from 0.0 to 17.5 percent, and on a No. 325 sieve the 
residue ranged from 0.0 to 24.2 percent. The fine- 


ness, as determined by the air-permeability appara- 


tus, ranged from 5,470 to 9,730 cm*/g. There was 
no relation between fineness and other properties of 
the cements. 

4. All the cements attained final set within 24 hr. 

5. The compressive strength of 2-in. cubes pro- 
portioned 1 to 3 by weight ranged from 310 to 1,980 
psi at 7 days, and from 690 to 3,160 psi at 28 days. 
The 12 masonry cements of portland cement-base 
type had strengths ranging from 1,000 to 1,980 psi at 
7 days and from 1,280 to 2,870 psi at 28 days. 
Strengths of mortars proportioned 1 to 3 by weight 
were slightly greater than those proportioned by 
volume. The 1-year tests indicated an appreciable 
increase in strength for most cements with both air 
and water storage. 

6. The water-repellency values for the mortars 
proportioned by weight ranged from 0.26 to 0.70, 
and from 0.32 to 0.76 for mortars proportioned by 
volume. The requirement of Federal Specification 


SS-C-181b for water repellency is that the ratio o 
the 1- to 24-hr absorption shall not exceed 0.90. 

7. Of the mortars proportioned by weight, only 
one had a water retention of less than 65 percent, 
and of mortars proportioned by volume, 4 had 
water-retention values less than 65 percent. 

8. The air content of 1 to 4 mortars ranged from 
3.7 to 29.0 percent. 

9. Three of the cements having less than 10 per- 
cent of air in the 1 to 4 mortar had a very low dura- 
bility factor. Only 2 of the 11 cements with more 
than 15 percent of air in the 1 to 4 mortar had a 
weight loss of more than 10 percent in 100 freezing 
and thawing cycles, and all had a very high dura- 
bility factor. 

10. Mortar bars of the different cements propor- 
tioned 1 to 3 by volume had shrinkage ranging from 
0.075 to 0.16 percent in 60 days. 

11. Neat pats made of these cements exhibited 
satisfactory soundness in both air and steam. The 
expansion of neat-cement bars as a result of auto- 


| clave treatment ranged from 0.03 to 6.2 percent. 


Fourteen of the cements had an expansion less than 
0.20 percent, the expansions of the other 3 were 0.56. 
0.64, and 6.2 percent. 
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Subsieve Particle-Size Measurement of Metal 
Powders by Air Elutriation 


Rolla E. Pollard 


Particle-size measurements of spherical metal powders by means of the Roller air an- 
alyzer, using samples up to 40 grams of powder, were found to be reproducible within approx- 


imately 
microns 
for particle diameters of 5 to 20 microns. 


tl percent of the original weight of the sample for particle diameters up to 120 
The separation limits closely approximated those required by Stokes law (V 


-=Kd 


For particle diameters of 20 to 70 microns, the 
size limits of fractions conformed approximately to the formula V 


:K''d'4, With particles 


larger than 70 microns, the accuracy of separation was uncertain. 
Reproducible results also were obtained with irregularly shaped particles, such as those 


of electrolytic copper powder. 


However, the accuracy of calculated size limits was uncertain, 


no attempt being made to check them by other methods of measurement. 


1. Introduction 


\ccurate and reproducible methods of determin- 

particle size are of great importance in powder 
metallurgy. The size distribution of the particles 
markedly influences the inherent packing behavior 
of a powder. This vitally affects the pressing and 
sintering techniques employed and the properties of 
the sintered object produced. 

Sieves are used almost universally for size analyses 

of particles large enough to be retained by a No. 325 
4443 microns). 
For smaller, or subsieve, particles, other methods of 
measurements must be used Among the many 
methods currently employed, air elutriation has a 
distinct advantage because dry, sized fractions can 
be ecllected and used for further study. 

One of the few commercially available devices for 

etermining particle-size distribution by air elutri- 
ation is the Reller air analyzer developed by P. S. 
Roller [1, 2, 3, 4]' at the U. S. Bureau of Mines. 
Tentative standard methods of particle-size analysis 
based on this apparatus are being considered by the 
Metal Powder Association [5], and the American 
Society for Testing Materials [6]. The Roller air 
analyzer has been used at the National Bureau of 
standards for several years in connection with studies 
ff the properties of metal powders. An investiga- 
tion was made to determine the effects of certain 
variations in equipment and operating technique on 
the reproducibility of the results obtained and the 
accuracy of the separations. These variations in- 

luded time of testing or criterion of end point, con- 

dition of interior surface of the settling tubes, par- 
ticle-size distribution of sample, range of particle 
sizes within a fraction, size of sample, humidity of 
air, and particle shape. Two prime objectives of 
the investigation were to establish the approximate 
range of application of Stokes law to metal-powder 
elutriation under the conditions found in the Roller 
analyzer, and to determine, if possible, the relation 
between the diameters and falling velocities for par- 
ticles of larger size. 


s1eVe 


cures in brackets indicate the literature references at the end of this paper* 
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2. Theoretical Background 


Particle-size measurement by the Roller air ana- 
lyzer is based on Stokes’ law, according to which the 
terminal velocity of a particle falling freely in a 
homogeneous fluid is given by the formula 


- (p—po)d’ 
Kg p p ’ 


where 


\'=terminal velocity, cm/sec 

K =a constant that depends on the shape of the 
particle and the nature of the air flow 

g=gravitational constant =980 cm/sec? 

p=density of the particle, g/cm* 

po=density of the fluid=0.0012 g/cm’ for air at 
21° C and 760 mm 

d=diameter of particle, cm 

u=viscosity of the fluid 
aic. 


1.82<10-4 poise for 


Stokes’ law is strictly 
cles and streamline 


valid only for spherical parti- 
flow. When A=1/18, then 
V=g(p—po)d*/18u or, with air as the fluid and d 
expressed in microns, V=0.00299(p—po)d*. Under 
nearly ideal conditions (that is, a single sphere falling 
in a relatively large cylinder of quiescent fluid) the 
experimental measurement of the rate of fall of a 
sphere is a convenient method of determining the 
value of the viscosity of a fluid. Castleman [7] has 
analyzed a number of such measurements to deter- 
mine the range of Reynolds numbers’ over which 
Stokes’ law gives accurate results. Castleman’s 
mean curve coincides with the curve representing 
Stokes’ law at approximately R=0.5. For iron and 
copper powders in air, this would correspond to 
diameters of 32 and 31 microns, respectively. 
However, such calculations assume ideal conditions 
not found in an air separator, such as the Roller, 
where large numbers of particles enter the settling 
? The Reynolds number, R, represents the ratio of the inertia terms to the 


frictional terms in the equations of motion of a fluid. In the case of a particle 
falling in a fluid, R=dVpo/u 





chambers simultaneously and the movement of 
individual particles may be influenced or interrupted 
by the presence of other particles or by the containing 
wall. The air also possesses motion, which changes 
in direction and in velocity while in contact with the 
particles and with the walls of the chamber More- 
over, most particles produced commercially are not 
true spheroids but are irregular in shape, thus fur- 
ther complicating such In this respect 
Wadell [8] has found that with nonspherical particles 
consistent results can be obtained im sedimentation 
or ant elutriation tests only whea the Reynolds num- 
ber does not exceed 0.2 (R<0.2 


tests 


3. Method and Materials 


The Roller air analyzer consists, essentially, of a 
device by which powder particles are blown into a 
settling chamber through which dried and metered 
air moves at a uniform measurable velocity. Means 
are provided for collecting the particles carried out 
of the settling chamber and for removing aad re- 
circulating lodged particles. A diagram of the appa- 
ratus used ia this investigation is shown in figure 1. 
Low-pressure air is cleaned in the dust trap, H, and 
the oil trap, G, before being dried in an activated 
alumina dryer, F. E is a vessel coataining a mois- 
ture-indicating compound. Air flow is regulated by 
means of calibrated metering orifices, A, and ma- 
nometers, Cand D. The air stream passing through 
the jet, J, at high velocity impinges on the powder 
sample, K. The U tube, M, holding the sample is 
given a jecky, oscillatory motion by bumper arms 
rotating on the shaft, L. The oscillating motion of 
the U tube (288 oscillations a min) continually 
brings fresh material under the action of the jet and 
newly exposed powder particles of amenable sizes 
are blown into the settling tube, P. Theoretically, 
particles smaller than the selected size are carried | 
out of the tube and collected in extraction thimble Q 


¥ 


\ 
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\ 
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FIGURE Details of Roller air analyzer assembly 
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During the investigation it was found that alu 
extraction thimbles, RAQS (coarse), were def 
superior to the paper thimbles originally used 
example, a paper thimble, during a 120-mi: 
with dried air” but without powders, will ga 
much as 0.81 ¢ in weight. This would amou 
over 2 percent of the weight of a 10 -g sa 
During this time an alundum thimble will gain 
0.004 g. Moreover the alundum thimbles ca 
cleaned readily in nitric acid or aqua regia and rv 

Means are supplied for dislodging and recirculating 
material that sticks to the walls of the settling 
chamber by a tapping device, N, with two leather. 
faced hammers that strike the tube about 870 timex 
amin. Later this device was supplemented 
Svyntron vibrator, O, (type V—4, 7,200 vibrations , 
min 

Four settling tubes with inside diameters of 9 
$'5, 24, and 1'4in., respectively, were supplied with the § 
apparatus. The three largest tubes vary in diamet, 
by the ratio of 1/2. Therefore, to the limit of appl 
cation of Stokes’ law (!"= Ad?) successive separations § 
made in these tubes at the same rate of air flo 
would yield fractions whose maximum parti 
diameter would vary in the same order by the rati 
of 2. The same rate of air flow used in the 1\%-i) 
tube would give ar air velocity lower than tha 
calculated by Stokes’ law by the ratio of (1\)?/(1 
Although it is not specifically mentioned by 
manutacturer, it is possible that this change in rat 
was intended as a correction for Stokes’ law at hig! 
air velocities 

The metal powders used in the investigation a 
listed in table 1. The permalloy powders, PC an 
PF, were found to be predominantly spherical | 
shape, 90 to 95 percent of the particles having near! 
perfectly circular projected images, as shown 
figure 2. A mixture of the two powders (2 parts b 
weight of PC to 1 part of PF, designated in table 
as PM) covered the subsieve particle-size rang 
fairly well and was used extensively in checking tly 
accuracy of particle-size analyses by microscopi 
measurements. The atomized copper and = zin 
powders were also largely spherical in shape 
shown in figures 3 and 4, respectively. approx 
mately 80 percent of the projected images of thes: 
two powders were either circular or slightly elliptica 
with ratios of minor to major axes of 0.9 or mor 
It was observed that there was not much differen 
in size distribution among the several fractions 
whether the perfectly circular images were measure 
or the minor axes of all elliptical images were i 
cluded. The controlling dimension in the separatior 
of these particles, therefore, seems to be the mino: 
aXiS 

The other powders of table 1 were more or less w- 
regular in shape. For example, the irregular shap 
of the electrolytic copper, Cu B, powder is shown 1 
figure 5. No attempt was made to check the par- 


) 


The term “dried air” refers, in this text, to air dried in an activated alu 


dryer to a nominal dewpoint range of —40° to —70 
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ticle-size dist ibution of such powders by microscopic 
measurements. They were used, however, in studies 
involving reproducibility. 
Attempts were made 

particle density of some 
according to ASTM D153-—39, “Standard Method of 
Test for Specific Gravity of Pigments”. The results 
were not very satisfactory, however, because the 
density of the kerosine changed during vacuum treat- 
ment to remove air occluded in the powder sample 
during its introduction into the kerosine. This 
change in density caused errors in the determination 
of the true volume of the sample. The range for 
permalloy powders, for example, was 8.70 +0.06 
gicm®, This uncertainty undoubtedly could be re- 
duced by employing suitable liquids of constant 
boiling point and possibly by using larger samples. 
It is believed, however, that the true solid density of 
all the powder particles used in this investigation, 
with the probable exception of the reduced tron 
powders Fe B, Fe B1, and Fe B2, are very close to 
the theoretical solid densities given in table 1. 


determine the solid 
the powders used, 


to 


4. Effects of Variations in Equipment and 
Testing Procedure on the Reproducibility 
of Results 


4.1. Time of testing—criterion of end point 

Separation by air elutriation, like separation by 
sieving, is never complete. One of the first problems 
encountered in using the Koller air analyzer, there- 
fore, is the adoption of some convention for estab- 
lishing an artificial end point. The method origi- 
nally used by Roller [1, 3, 4] was to operate by 30- 
min intervals, weighing the amount collected for 
each interval. His stopping point, or artificial end 
point, was reached when the weight of powder col- 
lected in any interval was less than 10 percent of that 
collected in the first interval. This method also is 
recommended in a booklet published by the Metals 
Disintegrating Co. [9], and is being considered in the 
tentative standards of the MPA [5] and the ASTM 
6]. The manufacturer, in the instructions sent with 
the apparatus, recommends a similar procedure but 
with 2-min intervals. 
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Later Roller [10], working with some ceramic ma- 
terials, found that when the weight of fraction sepa 
rated was plotted against the rate of separation, » 
straight line or characteristic was obtained, tha: 
when extrapolated to zero rate, gave the total amount 
of material of the selected-size range originally pres- 
ent in the sample. The time intervals used fo 
portland cement ranged from ' to 25 min. Rolle 
[11] has recently applied this end-point method t 
the particle-size analysis of metal powders, using 
time intervals as short as 1 or 2 min, which, in th 
light of the results of the present investigation 
appear to be of too short duration. 

For example, a number of tests were made with 
the permalloy powder mix, PM, using time intervals 
of 2, 4, 8, 15, and 30 min. Some of the tests wer 
stopped after the interval in which the amount col- 
lected was less than 10 percent of that collected i 
the first interval (10-percent end point). Others 
were continued until the collection per interval was 
less than 1 percent of the first interval (1 percent end 
point). 

In figure 6 the amounts of the 0- to 10-micror 
fraction collected during some of these tests ar 
plotted against the rate of separation, according to 
Roller’s method [10, 11]. It will be noted that thos 
portions of the curves representing the early stages 
of the test (which were used by Roller) tend to 
flatten out with decreasing time intervals. In the 
2-min-interval curve, the 4 points representing the 
third, fourth, fifth, and sixth intervals may be fairly 
well represented by a straight line. Extrapolation 
of this line to zero rate, however, indicates a rather 
low result, 9.55 g, or 21.4 percent, of a 40-¢ sample 
as compared to 27 +1 percent as measured by the 
10-percent end point. 

This may be due, at least in part, to the fact that 
for very short time intervals the largest particles of 
the selected size range do not have time to leave the 
settling chamber. The resultant velocity, Vx, of a 
particle in a settling tube depends in magnitude and 
direction on the difference between the velocity of 
the air, V,, and the free-falling velocity of the par- 
ticle, Vr. The 4%-in. settling tube, for example, 
has a section of uniform diameter about 28 cm in 
length. In order for a particle to escape from this 
tube Vy must be greater than 28/60¢ cm/sec, where 
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rmalloy powder (p=8.7) with an air velocity 
lated for a nominal particle diameter of 10 
ons, it can be shown that the largest particle 
that could escape from the tube in a given time 
val would be 9.1 microns for 1 min, 9.5 for 2 
min, and 9.97 for 30 min. For very short time 
tervals, therefore, the actual separation would be 
at a diameter appreciably smaller than that for 
vhich the air velocity was calculated, and the results 
would be low, particularly if the size distribution of 
© the sample were such that large numbers of border- 
line particles would be affected. 

The influence of the time interval on the end- 
point criterion is illustrated in figure 7, which repre- 
sents all of the tests made on the 0- to 10-micron 
fraction of the PM powder. The upper and lower 
solid curves represent, respectively, the highest and 
owest values obtained with any time interval 
within the total testing time indicated on the scale 
of abscissas. The intermediate broken curves rep- 
resent the lowest values obtained with 2-, 4-, 8- and 
|5-min intervals. The total number of tests include: 
2-, 4-, and 8-min intervals—4 tests each, with 1 in 
each case carried to the 1-percent end point; 15-min 
ntervals—S8 tests, with 3- to the l-percent end point; 
‘0-min intervals—67 tests, with 3- to the 1-percent 
end point. 

For tests with time intervals shorter than 30 min, 
test conditions were uniform, except that in those 
continued to the 1-percent end point a Syntron 
vibrator, O, figure 1, was used in addition to the 
regular hammer, N, figure 1. The 30-min-interval 
tests, however, included many in which the test con- 
ditions were varied considerably. For example, the 
three points just outside the “‘knee’’ of the upper 
curve (at 30, 60, and 90 min) were obtained without | 
a vibrator but with humidified instead of dried air. 
Other tests were made with and without the vibrator 
and with samples varying in size from 20 to 100 g. 
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It will be noted that, for the 30-min interval, in 
spite of a very wide variation in initial separation 
rate (from 13 to 26 percent per interval), the results 
tended to converge so that the 10-percent end points 
(marked with barred arrows) varied no more than 
+1 percent of the original weight of the sample. 
For the shorter time intervals the 10-percent end 
points were progressively lower. These results show 
conclusively that to get within the end-point range 
of the 30-min intervals, tests with shorter intervals 
would have to be carried beyond the 10-percent end 
point. All the 1-percent end points (marked with 
plain arrows) were within this range except the 2- 
min interval, which was at the lower edge. Roller’s 
end point (as determined in fig. 6) would be at 21.4 
percent on the upper curve of figure 7. 

It was concluded from these results that 30-min 
intervals and the 10-percent end-point criterion 
were most satisfactory. This procedure, therefore, 
was adopted as standard for all subsequent tests. 
Complete Roller analyses of permalloy and other 
metal powders are given in tables 2 to 9 inclusive. 
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TABLE 4 Roller analyses of hydrogen reduced tron powders, Fe Bl and Fe B2, 40-¢ samples, minus No. 325 sieve 30-min 
10-percent end point, dry air 
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4.2. Condition of Interior Surface of Settling Tube 


The tabulated results show that for diameters up 
to 20 microns, good reproducibility (on the order of 
+ |] percent of the original weight of the sample) 
was obtained with powders such as Cu A, Cu B, Fe A, 
Fe B1, and Fe B2, which contained relatively small 
amounts (5 percent or less) of the 0- to 5-micron 
fraction (tables 3. 4. 5 Good results were also ob- 
tained with the PM powder (0- to 5-micron fraction, 
theoretically, about 13.5 percent) without the use of 
a vibrator if the initial separation was made at 10 
microns (series 3, table 2). However, for separations 
made at 30 or 40 microns in the 1's- and 1 '-in. settling 
tubes as received, the results often varied widely. 
These variations were attributed to the rough con- 
dition of the interior surface of some of these tubes. 
The 9-, 44-, and 2\-in. tubes were made of welded 
stainless steel sheet, and the interior surfaces were 
smooth. The original 1-in. tubes supplied with the 
unit were made of brass tubing and failed after 
relatively short periods of service. One of these 
tubes was used only in series | of table 4, and gave 
reproducible results; another was used only for 
series | of table 5 and gave variable results. Whether 
or not these inconsistencies were due to the condition 
of the interior surfaces of these tubes cannot be 
stated, because they were not examined upon re- 
moval prior to returning them for replacement. 

In all the other tabulated series of tests requiring 
the use of small tubes, 1'- or 1's-in. diameter, aus- 
tenitic stainless steel tubes were employed. Both of 
these tubes were made of stainless-steel tubing, 
machined to the required inside diameter, and the 
interior surfaces were relatively rough. Powder 


particles adhering to these rough surfaces were 
difficult of 
undersized material were often carried over into the 
next 


microscopic measurements of the 30- to 40-micron 


to dislodge, and substantial amounts 


fraction. In series 3 of table 2, for example, 


fraction showed the presence of 24.8 percent of | 


undersize material, which had obviously remained 
as a residue from the prior 30 micron separation. 
Such condition caused wide variations in the results 
obtained for both the 20- to 30- and 30- to 40-micron 
fractions. The of 
was considerably reduced with consequent improve- 


carryover undersized material 
ment in reproducibility when a vibrator was used in 
addition to the regular hammer (series 4, table 2). 
Still further improvement in this respect was ob- 
served when the 30-micron separation was made in 
the smooth 2\-in. tube 
the interior surfaces of the smaller tubes were lapped 
2). The amount of oversized material 


(series 5, table 2), or when 
(series 6, table 
was reduced materially only by lowering the air 
velocity (series 7, table 2). Reproducible results 
were obtained with all other powders tested with 
these tubes when the interior surfaces were lapped 
(tables 3, 4, 5, 7, and 8). 


4.3. Dislodgment of Small Particles 


Efficient separations at the 
powders such as PF, PM, Fe C, Cu C, an 
which contain substantial amounts (13.5 to 
percent) of 0- to 5-micron fraction, were fou 
be difficult because of the tendency of partic 
adhere to the walls of the settling tube. In 
cases the hammer device alone was not suffici 
effective in dislodging the particles, consequent! 
results were variable and the total recovery wa 
(series 1, tables 2 and 6). But when a vibrato 
used in addition to the regular hammer device, \ 
tions in results were diminished and the tota 
covery was markedly improved (tables 2, 6, and 7 


5-micron ley 


4.4. Size Range of Fraction 


It will be noted from the tables that a separatio 
at a larger particle diameter usually gave somewha 
lower results than the accumulated weights of th 
finer fractions. In series 2 of tables 2 and 3, fo 
example, the average for the 0- to 20-micron frac. § 
tions of PM and Cu B powders was 68.9 and 49) 
percent, respectively, whereas the sums of the fin 
fractions (0- to 10- and 10- to 20-micron in series 
table 2, and 0 to 5, 5 to 10, and 10 to 20 in series | 
table 3) were 70.8 and 53.4, respectively. This 
believed to be a consequence of the uncertainty « 
the end point under some conditions. For the 0- t 
20-micron fractions, the rate of separation for th 
initial 30-min intervals was very high (as much as 
60 percent by weight of sample). However, this 
rate fell off rapidly, and at the 10-percent end poin 
as much as 4 or 5 percent was still being blown ov 
per interval. This clearly indicates that some o 
the border-line particles were not being separated 

Similar end-point uncertainty was found wher 
accumulated, sized, fractions were analyzed by th 
Roller apparatus, as shown in table 9. The samples 
in these tests were the accumulated fractions of th 
5 tests included in series 7 of table 2. These frae- 
tions contain many border-line particles that ma 
be blown over or remain in the settling tubes, de- 
pending on the number of such particles present an 
the rate of separation. A sample of accumulate 
20- to 30-micron fractions, for example, may con- 
tain some 18 and 19 microns as well as 31, 32, an 
possibly a few 33-micron particles. If a separatior 
is made at 20 microns, the initial rate of separatiot 
will be low but will fall off slowly. so that by th 
time the end point is reached all of the undersiz 
most of the 20-micron, and many 21- and 22-microt 
particles will have been blown over. This may 
result in the separation of a nominal minus 20-micro! 
fraction of appreciable size, even though repeated 
microscopic measurements have indicated that 9 
to 99 percent by weight of the original sample was 
within the 20- to 30-micron size limits. When th 
separation is made at 30 microns, however, the initial 
rate of separation will be very high and will fall off 
rapidly, so that some of the border-line particles 
will remain in the tube to be blown over during th 
40-micron separation. 
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Sand with hammer only. 


4.5. Size of Samples 


size of sample, within the relatively small lim- 
posed by the diameter of the sample tube (M, 

appeared to have little effect on Roller par- 
ze analyses. This has been indicated pre- 

by the inclusion of tests (3 for each weight) 
ples, weighing 20, 30, 40, 60, 80, and 100 g, in 
sults represented by figure 6. It was noted, 
er, that the results for the 80- and 100-g samples 
ear the lower limit of the range, and large num- 
f such tests might broaden the range. On the 
hand, with very small samples, particularly of 
»wders, the small losses of material unavoidably 
ed become increasingly important on a per- 
ve basis (see total recovery values of tables 2, 
nd 8). In this investigation, consistent results 
obtained by following the suggestion of the 
ial Disintegrating Co. [9] in using samples weigh- 
bout 10 times the numerical value of the appar- 
density of the powder, with a maximum of 40 g 


iost cases. 


t 


4.6. Humidity of Air 


lhe tendency of small particles to stick to each 
other or to the walls of the settling tubes was found 
to be markedly affected, in some cases, by the humid- 
ty of the air supplied to the settling tubes. The 
solid curves, E, A, and B shown in figure 8 represent 
the range of results obtained during the separation 
of the 0- to 10- or 0- to 5-micron fractions of the PM, 
Fe C, and Cu B powders with dried air and with 
hammer only. The broken curves E’, A’, and B’ 
represent the results obtained with humidified air 

The air was humidified by 

passing it, in two stages, through alundum thimbles 
mmersed in large bottles of water. The relative 
humidity, as measured by wet- and dry-bulb ther- 
mometry, was about 90 percent 

Curves E and E’ of figure 8 show that high humid- 
ty greatly increased the separation rate of the 0- to 
\0-micron fraction of the PM powder. The final 
results, however, were within the range of those ob- 
tained with dried air. The effect of humidity on the 
initial rate of separation of the 0- to 5-micron fraction 
of the Cu B was negligible, but the final results were 
considerably higher than those obtained with dried 
air (curves B and B’). Humidity had little effect on 
carbonyl iron powder, the results at all stages being 
within the range of those obtained with dried air 
curves Aand A’). With this powder, and with Cu C 
powder, more efficient dislodgment of particles was 
obtained by the use of a vibrator in addition to the 
hammer as shown by curves C and D (see also tables 

6,7, and 8). 

lt is probable that the varying effect of high hu- 
midity may be connected with the static charges 
picked up by the particles during separation. These 
charges may vary in nature and magnitude with the 
humidity of the air, as well as with the size and com- 
position of the particle. 


AMOUNT COLLECTED,% OF ORIGINAL SAMPLE 


60 36 80 300 
30 é 80 270390 
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Ficure 8&8 Effect of the humidity of the air stream on the 
separation of jine fi actions of metal powders 


Curve E, 0- to 10-micron fraction of permalloy-mix powder obtained with dried 
iir without vibrator; E’, same conditions except for relative humidity of the air 
about 90 percent; A, 0- to 5-micron fraction of carbonyl iron, Fe C, powder ob- 
tained with dried air without vibrator; A’, same except for relative humidity 
of the air—about 90 percent; B, 0- to 5-micron fraction of electrolytic copper, 
Cu B, powder obtained with dried air without vibrator; B’, same except for rela- 
tive humidity of the air—about 90 percent; C, 0- to 5-micron fraction of carbonyl! 
iron, Fe C, powder obtained with dry air and with vibrator and hammer 
1), 0- to 5-micron fraction of electrolytic copper, Cu C, powder obtained with 
dry air and with a vibrator in addition to the regular hammer 


4.7. Particle Shape 


With the metal powders used in these experiments, 
particle shape was found to have no effect on the 
Roller analyses insofar as reproducibility was con- 
cerned. With the modifications in equipment pre- 
viously discussed, and with reasonably uniform 
operating conditions, reproducible results were ob- 
tained with all powders tested at air velocities cal- 
culated for diameters of 5 to 40 microns. However, 
the accuracy of these separations, or rather the 
accuracy of the values assigned to the size limits 
of the fractions, depends on the relation that is 
assumed to exist between the diameter of the particle 
and its falling velocity. The accuracy of some of 
these assumptions will be discussed in the next 
section. 


5. Accuracy of Separations 


The accuracy of Roller analyses was checked by 
microscopic measurements on sized fractions of perm- 
alloy, atomized copper, and atomized zine powders. 
Samples of these powders, dispersed in xvlene con- 
taining a few dreps of cedar oil on plain glass slides, 
were mounted on a microscope; the images were 
projected with transmitted light upon a screen at 
1,000 magnifications; and measured to the nearest 
millimeter. While it is recognized that there are 
inherent errors in this method of measurement, they 
are minute with respect to the present application 
of the method and, therefore, may be overlooked. 

Microscopic measurements indicated that sepa- 
rations of these nearly spherical particles closely 
approximated the requirements of Stokes law for 
diameters of 5 to 20 microns. The parabolic curve A 
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of figure 9. for example, shows the theoretical values 
for the terminal falling velocities of spherical PM 
particles, as calculated by Stokes law, plotted against 
the particle diameter in microns. In actual separa- 
tions at air velocities taken from this curve, it was 
found that a minimum of 98.5 percent by weight of 
the 0-to 10-and 10-to 20-micron fractions were within 
the designated size limits, as shown by the micro 
Similar results 
were obtained with 0- to 5-, 0- to 10-, 0- to 20-, 5- 
to 10-, and 10- to 20-micron fractions of Cu A, PF, 
Zn, and PC powders, as shown in tables 3, 7, and 8. 

For separations at larger diameters, air velocities 
calculated by Stokes law were too high. For ex- 
ample, the air velocities required for separations of 
PM powder at 30 and 40 microns would be 23.4 and 
41.6 «m/sec, respectively (see curve A, fig. 9). It 
was found, mostly by trial and error, that accurate 
separations of this powder could be obtained with 
air velocities of 18.7 cm/sec for 30 microns, and either 
27 or 28 cm/sec for 40 microns. 

Dallevalle [12], citing the work of Allen [13] and 
Martin [14], states that, between the region of 
streamline motion, where V varies as the square of 
the particle diameter, and the region of fully turbu- 
lent motion, where V varies as the square root of 
the particle diameter, there is a rather extensive 
region where a straight line relationship exists and 
V varies as the first power of the particle diameter. 

Curve B of figure 9 is a straight line through the 
points 18.7 and 27 determined experimentally as 
suitable for accurate separations at 30 and 40 
microns respectively. This line cuts the parabola 
of Stokes law at about 20 microns. Tests were 
made of PC powder, using air velocities taken from 


scopic counts of series 3 of table a 





30 


curve B for separations at diameters rangi: 
10-u increments from 40 to 90 microns (40. 
ete, These velocities were found to be too lo 
accurate separations. The nominal 70- to 80-n 
fraction, for example, when measured micro 
cally actually was 100 percent within the 6 
70-micron range. Microscopic counts were ma 
each of the fractions to determine approximate! 
actual separation points. These points, when pk 

on logarithmic coordinates against the air velocities 
used, fell near a straight line with a slope of approx 
mately 1.4. 

This suggests a relationship similar to that 
by Burke and Plummer [15], who related their ow 
and other available data to the general formula fo 
the resistance to motioa, R, of a particle in a flyj 
R kh po" 12" (d/2)" The force causing Motor 
is for a sphere, F'=(md*/6) (p g. At termina 
velocity R= F and 


Do 


(2) 


"rg 


6k ‘« 


y" (p— Po) po’ 


When l, 
streamline motion, and when n 
Newton’s law for turbulent motion. 

Burke and Plummer suggest that in the inte: 
mediate region I< n<2, and probably changes 
smoothly and continuously betweea the two limiting 
values. 

If n is made equal to 1.25, eq | 


[ ] Po 0 uu 0 7 he + 


With k= 16.07, this gives for permalloy powde 
V==18.7 at 30 microns and V=27.9 at 40 microns 
A theoretical curve for PM powder with air velocities 
calculated by this formula is shown as curve C i 
figure 9. 

Tests were made on unsieved PC and Cu A pow- 
ders with air velocities calculated by eq (3) fo 
particle diameters larger than 20 microns. The re- 
sults shown in table 8 were comparable to micro- 
scopic measurements for particle diameters up to 70 
microns, & minimum of 94.6 percent by weight o! 
the 60- to 70-micron fraction being within the spec- 
ified size limits. These results indicate that, for 
particles larger than those which obey Stokes law 
the exponent n of eq (2) changes very rapidly from 
a value of | at 20 microns to about 1.25 at 30 microns 
Between 30 and 70 microns the change in n, if any 
is too small to be detected by this method of testing 
Above 70 microns the separations became progres- 
sively less sharp. There appeared to be some pre- 
ponderance of oversized material in the larger frac- 
tions that may indicate that n changes in this region 
to a value larger than 1.25. However, an increasing 
overlapping of size range among consecutive frac- 
tions (both undersize and oversize) made it impossibl: 
to draw definite conclusions. In spite of this over- 
lapping, the results continued to be reproducible up 
to 120 microns, as may be seen in series 4 of table § 
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previously stated, no attempt was made to 
the Roller analy Ses of powders composed of 
arly shaped particles by microscopic measure- 
Some work has been done in this field, but 
present time, the exact nature of the dimen- 
hat control the measurement of an irregularly 
particle by air elutriation is not well under- 
or is the possible relationship of these dimen- 
pany v isible dimensions or to the surface area 
ne of the particle 


6. Summary 


eve particle size analyses of a number of 
sowders were made with a Roller air analyzer, 
(-min. intervals and a 10-percent end-point 
mn This procedure was found to be 
for samples of normal size distribution. The 
int was somewhat uncertain, however, when 
initial rates were very low but persistent, as 
paration near the lower limit of an accululated, 
fraction, and (2) the initial rates were very 
ut fell off rapidly, as in a separation near the 
limit of an accumulated sized fraction, or with 
itively large range of sizes within the fraction 
as 0 to 30 microns 
Particle shape appeared to have no effect on the 
reproducibility of the tests Results reproducible 
within approximately +1 percent by weight of 
sample were obtained with very irregularly shaped 
particles for diameters up to 40 microns, and with 
rical particles for diameters up to 120 microns 
provided the following conditions were fulfilled: 
{n alundum extraction thimble was used to 
ollect the fractions. 
2». The inside surfaces of the settling tubes were 
ooth 
\ Syntron vibrator was used in addition to the 
cular hammer, particularly for powders containing 
Jatively large amounts of the 0- to 5-micron 
raction 
t. The size of the sample, while not cirtical, was 
maintained at approximately 10 times the apparent 
density as a lower limit, with 40 ¢ as a maximum in 


satis- 


most cases 

The effect of high humidity on the separation of 
metal powders was not consistent. With some 
powders the rate of separation was markedly in- 
but with others no effect was noted. In 
view of this uncertainty, it would seem advisable to 
use dried air or air of relatively constant humidity 
for all tests 

The accuracy of analyses of spherical powders was 
checked by microscopic measurements of the frac- 


creased, 


tions. With these powders Stokes law applied for 
diameters of 5 to 20 microns, and the air velocity 
could be calculated accurately as proportional to the 
square of the particle diameter (V=Kad*). For larger 
diameters, the air velocity was found to be propor- 
tional to the 1.4 power of the diameter (V=A’'d'' 
At still larger diameters, increased overlapping of 
size range among consecutive fractions prevented 
any definite conclusions, except that the results 
continued to be reproducible 


Grateful acknowledgment is made to W. Earl 
Lindlief, formerly of the National Bureau of Stand- 
ards, who assisted in the planning and supervision 
of the early part of the investigation; and to John C 
Everts and Caroline Plummer Jeffries for their 
patient and efficient assistance in the monotonous 
and time-consuming tasks of making the vitally 
important particle-size counts and analyses. 
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Eigenvectors of Matric Polynomials 
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2 
4 
2 
ie It is the object of this paper to compare the eigenvectors of an arbitrary nn matrix A 
4 over the complex field with those of the matrie polynomia! f(A While it is well known 
that each eigenvector of A is an eigenvector of f(A), it is not, in general, true that A and 
f(A) have identical eigenvectors. In this regard a necessary and sufficient condition that 
A and f(A) have identical eigenvectors is given. The condition is that both (1) and (2) hold 
l f’(x) #0 for all eigenvalues \ of the matrix 4 corresponding to nonlinear clementary 
livisors 
2) The values of f(z) are distinet for all eigenvalues w of the matrix A corresponding 
to linear elementary divisors 
When either (1) or (2) fails to hold, then f(A) has eigenvectors that are not eigenvectors of 
1. This situation is also discussed 
vector space of eigenvectors of the matrix A | then 
esponding to the eigenvalue nN shall be denoted _ " 7 
Aj. The vector space spanned by the eigen- ) 3! m—1)! 
tors of A shall be denoted by [A]. Further, 
A and d } [A] shall denote their dimensions 0 . ‘ 
. . . ° . >" )\1 
lt clear that each eigenvector of A is an eiwen- = ™ as 
f vr of f(A That is, 1,[A] = Von (7 A )] for each 
nvalue XNof A. Thus VLA] V[fCA)] 0 0 3 
wt J be the Jordan canonical form of A. Then | ¢¢7 
Oe thy exists a nonsingular matrix P such that 
BP AP=J and so P-'(A)P=fJ 
Lemma | The C1 ue nvectors of A and f(A) are iden- 
f and only tf. the eige nvrectors of J and f(J) are 
ul Proof. This follows from the fact that PE is an 0 0 ’ ’ 
4 envector of PLP if € is an eigenvector of the | oO 0 0 0 0 . J 
B matrix i 
Since Pi VAI I s thes “— : , 
| ‘ l ; J A), a follow th a Proofs of the lemma miay be found inh Wedderburn 
; : 


V{f(A)]} —d{ VLA}} =d{ VIP} —d i VI), 
here Pi V\It denotes the space of all vectors P, 
is an eigenvector of J corresponding to \ 





hey é 
re ¢& 


: Lemma 2? If dD diag la.a,. » a en so ae 
¢ T.%, tr], where a.B. nw are distinet. then 
: D Vf D))) Furthermore, the ¢ igenvectors of D are 
. tical with those of {(D if, and only if. f(a), F(B). 
f(r) are all distinet 

Proof. Uf Dis of order n, it is clear that both the 


venvectors of D and f(D) each generate the whole 
dimensional vector space. If f(@)=/(8), then it is 
asily seen that /(27) has eigenvectors corresponding 
to the eigenvalue f(a@)=f/(8), which are,not eigen- 
ctors of D. 
Lemma 3 If 


a es @ 0 0 0 


0 d l 0 UO 0 


SO Peas BPA SE 


i109 00..- DO QD a 
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and MacDuffee 


Lemma 4 If f’(a)=—f"" (a) f?(’)—0 and 

/ h) +0 (7 i. ee m—2). then 
df Veaol f(Im(d))]} —df Valen (A) 
and conversely Also. 
di Vea flIm(X))] Li Vy Imn(X)] =—m—1 

if, and only if, POX / Ay=O0 

Proof If f’(r f’’(r f(r 0, but 
f"* (x) <0, it follows that the rank of the matrix 
f(J»())—fONL is m r+1), as may be observed by 
noting the nonzero minor of order m—(r-+-1) in the 
upper right-hand corner of the matrix. Thus the 


nullity of f(J,0))—fOOL is r+1. Tf the first m-1 
derivatives of f(z) vanish at a=), then the nullity 
of f(T m(X))— FONT is m. In a similar fashion one 
finds that the matrix J,(4)—A/ has nullity equal 


to 1. This establishes the lemma 

Corollary. The matrices J,,(4) and f(.J,,(4)) have 
identical eigenvectors if, and only if, f(A) #0 

It is to be noted that in this case Vea [f(J,())] 


23. H. M. Wedderburn, Lectures 


on Matrices, (Am. Math. Soc., ¢ 
1ui4 


Noquium 


Pub. vol. 17, New York, N.Y 
Cc. C. MaecDuffee, An Introduction to Abstract Algebr John Wiley & Sons, 
New York, N. Y., 1940). 


















and | {J 
vector [1,0, 0] 

At this point it will be convenient to introduce the 
language of elementary divisors. Each block in the 
Jordan canonical form of a matrix A corresponds to 
an elementary divisor (r—X)" and conversely. 
Each such block is called the hypercompanion matrix 
J,() of the polynomial (2— A)” 

Suppose the nonlinear c 


\)} are each spanned by the single column 


T 


of 


divisors 





ementary 


the matrix A corresponding to the cigenvalue 
are of the form (*4—\)”™ r—r)™ , , (r—A)™*, 
where m; i=1,2. ky) are integers such that 
mi; mi: mi Furthermore, suppose 


that («—A appears as an clementary divisor n; 


a" 
l 
nonlinear elementary divisors of A corresponding to 
the eigenvalue X Denote the number of 
linear elementary divisors of A corresponding to the 
eigenvalue Further, denote by Ay the direct sum 
of the hypercompanion matrices of the elementary 


times Let —Pr denote the total number of 


by 7 


divisors corresponding to A. Set 
dy=di Vea, [f(A] ) —d) [Ag] 
Lemma 5 Tf the first m 1 derivatives of f(x) 
vanish at x=). then d, Son m | whereas, if 


l 
the first 7 mo l derivatives of f(r) vanish at a X. 


i 


but f°) (\) 40, then dy Sn} / > > * (ny, 1), 
i= | i 
where m" is the large st of the inte gers my, 2=1.2.....h,) 


such that m? <7 

Proof. This follows from lemma 4, the fact that 
the nullity of A,—X/ is the sum of the nullities of 
the characteristic matrices of the hypercompanion 
matrices of the individual elementary divisors corre- 
sponding to the eigenvalue \ and a similar statement 
about the nullity of f(A) —f(rayl 

Corollary. The matrices Ky and f(A) have identical 
¢ ie nvectors if, and only if, f’(\) #0 


Proof. This follows from the fact that d,—0 

The Jordan canonical form J=diag [K.,Kg,. . ..A,] 
of A is a direct sum of matrices Ay, where \=a.8, . . 
runs through the distinct eigenvalues of A. The 


subsequent theorems of this paper involve the fol- 
lowing main conditions: 

Condition l f’(r) Zé) for all ¢ rie nevalues of the 
matrix A to cleme ntary 
divisors 

Condition 2 The 
eige nvalues u of the matrix A corre sponding to linear 


corresponding nonlinear 
values f() are distinet for all 
elementary divisors. 

The first theorem concerns the case where condi- 
tion | does not hold, whether condition 2 holds or not 


If condition 1 does not hold. then 


d{V{[A] 


Theorem 1. 


di V{f(A)] } Sod 


\ 


Ay 


where X varies through all distinct eigenvalues of A 


corresponding to nonlinear ele me ntary divisors, and 
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d, is computed as in lemma 5. 
Proof One notes first that 





VJ] =ValJ] + Val) V5[J], 


and 





VAT = Vw (MD) + View lf] + Vee [ fd 
If f(a), f(8), . . ., f(a) are distinct, for a fixed « 
value A, the nullity of J—Al the same 
nullity of A,—AJ, and the nullity of f(J) 
the same as the nullity of f(A\)—fQ)l. It 
that 


is i> 


AOI 
fo 





dV [yal f(A)]} —di Vl] =ds 


Summing over distinct eigenvalues, one obtains , 


di VUf(J)]} —di VW] 


Dod. 
By the statement following the proof of lemma 
follows that 

di VUf(A)] 


d{V[A]} =o 


\ 


AY 


If 
f(J) 
J—sl, . 


f(a)=f(B ; f(p), then the nullity 
f(a)IT minus the sum of the nullities of J—a/ 
, and J—pl is d,+ds+...+d, and 


A Viel f(J)}}—df Veld]} —di Vs 





a. d{ VJ 


Summing, one obtains in either case 


d{ V[f(A)]} —d{ VA] 


D4, 


\ 


where \ runs through all the distinct eigenvalues « 
A corresponding to nonlinear elementary divisors 

The next theorem concerns the in whieh | 
condition | p 


- —-_ 


Cause 





holds but condition 2 does not 


Theorem2. If condition | holds,than V[A]=V{f(A 
In this case dd, 0, where X varies over all dis 
d 


tinct eigenvalues corresponding to nonlinear elemet ' 
tary divisors. 
If f(a) =f(B) f(p), 


di VJ] 


suppose that 
d{ Va [f(J)] 


t-d{ Vs[-J] 


J 


Ld{ V [J] 





Then as in lemma 2 the s-dimensional 
VulJ]+ Vl] +. . .+1,[J] contains vectors that ar 
not eigenvectors of J, but VV, - [f(J)} consists soli 
of eigenvectors of f(J) corresponding to f(a) 
Corollary. If f(a) 4(B) t(p), then 


spa 








df Vpa){ f(A)]} —d{ VLA] 


Ny), 


© 
4 
) 


Dod + (pr 
\ 








where runs over all eigenvalues among a, 8, 
wh ich corre spond to nonlinear eleme ntary divisors « 

























Dy n, is the total n umber of elementar y dirisors 
neal and linear) corre sponding tor 


allv. the next result covers the case in which 


tions 1 and 2 both hold 


orem 3 The matrices A and f(A) hare identical 
tors if. and only if (1) f(A) 40 for all ¢ hye n- 
\ of the matrix A corre sponding to nonlinear 
lary divisors. and 2 the values dis- 
y all eige nvralues u of the matrir A corre spond) iq 


fiw) are 


elementary divisors 
ark 1. It may readily be seen that VJ.J] and 
each can be generated by a set of linearly 
ndent vectors, each of which has a 1 in a single 





onent and 0 in all the remaining components 





J 
or 
| 
’ 
\ 
J 
; 
= 2 
Is 
his 
} 
dis. 
net ; 


12h \ 


par 








init 


if J—P AP is the Jordan canonical form of 
in VIFCA can be venerated by a subset of the 
vectors of P 
ark 2. A simple application of the foregoing 
shows that if A is a 22 matrix, then A and 
ive the same eigenvectors unless either l A 
onable and has distinct eigenvalues a, 8 for 
fla f(8), or (2) A is nondiagonable and 
] where a is the eigenvalue of A. 
an arbitrary polynomial, and ¢ ts an arbitrary 


’ 


L a) GC. 


following example is given as an illustration 








I l 0 0 0 0 0 oj 
0 0 l 0 0 0 0 
0 0 0 | 0 0 0 
{ 0 0 0 0 0 0 0 
0 0 0 0 0 0 0 
l l I l l 2 l 
= 0 i | 0 l | 
l ! } ! 2 








be verified that 


' l 0 0 0 0 0 0 
0 I 0 0 0 0 0 
0 0 l 0 0 0 0 
0 0 0 0 l 0 0 








and that the Jordan canonical form of A is 











ri Oo 0 0 0 0 oO 
00 1 0 0 0 0 
0 0 O 0 0 0 
j= P-i AP 00 0 0 0 0 0 
00 0 0 0 0 0 
0 0 0 0 0 | | 
Lo 0 0 0 0 0 1_| 
where 
/ 
ro 0 0 
| 01 0 
0 0 0 
A L.0 0 0 Lo t) 
00 0 0 
0 0 0 
00 0 0 
| l 
A / l 
= 0 l 
Since 
( 0 | | 2 
0 ] 0 1) <0) 


follows that 








It is readily that V[J] is generated by 
[1,0,0,0,0,0,0]", V,[-7] is generated by [0,1,0,0,0,0,0] 


seen 


and [0,0,0,0,1,0,0]7, and V_,[-J] is generated by 
[0,0,0,0,0,1,0]! Since VIJ]= Vi[J]+ VolZ}4- Vf] 
it follows that d/ [J] $. Since f|.J] is a scalar ma- 
trix,d{ VIf(J)]}}=7. Hence di{f(A)]} —d{ VIA} =3 
The same result is arrived at by the use of theorem 1, 
where S3d, is calculated as in lemma 5. Since 
» 
f'(0)=0, f’’(0)=0, and f’(—1)=0, where 0, —1 are 
the eigenvalues corresponding to nonlinear elemen- 
tary divisors, it follows that \ 
d Neo ilo l ]-2 2 


and 


From a glance at the above eigenvectors gener: 
\ [J], it is clear that VLA] is generated by the 
second, fifth, and sixth-column vectors of P. s 
F(A) is a sealar, V{f(A)] is generated by all s 


, , — - stidiinae 9 oe 
as before. By the corollary following theorem 2, column vectors of P- 


The author thanks A. H. Clifford, O. Taus 
and N. Wiegmann for their helpful suggestions. 








which checks with the observed results above WasHINGTON, December 1, 1952 
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The System Barium Oxide-Boric Oxide-Silica 


Ernest M. Levin and George M. Ugrinic 


A phase equilibrium diagram for the system BaO-B,O,-SiOQ, has been constructed from 
data, obtained essentially by the quenching method, on 178 ternary compositions. One 
new compound, 3Ba-O3B,O,-2S8i0,, melting at 1,009° C and possessing an extremely flat 
primary field appears, and its optical properties and X-ray diffraction data are given. The 
BaO-4B,0,; and SiO, primary phase areas include almost 70 percent of the diagram, and 
over these areas exists a large region of two-liquid immiscibilits Limited data indicates that 
the immiscibility gap decreases with increasing temperature. The solid region 
between 2BaO-3SiQ0, and BaO-2SiO, extends into the interior of the diagram and is complex 
in nature. An isofract diagram for the quenched glasses is shown. For precise tempera- 
ture control of a quenching furnace, a self-adjusting, bridge-tvype controller is described 





solution 


briefly 


l. Introduction 


(he system barium oxide-boric oxide-silica is of 
indamental importance in the ceramic field, and the 
three end members comprise the major constituents 
Sof the large class of barium crown glasses. Other 
jmportant constituents maybe alumimum- or zinc- 
xide in the dense barium crown glasses, lead in the 
barium flints, and alkalis in the light barium crowns 
v flints [1].!. Barium oxide surpasses nearly all 
ther oxides in imparting both high index and low 
lispersion to glass. Many of the well-corrected, 
present-day optical systems, such as are found in 
shotographic and projectional equipment, have re- 
ulted largely from the use of high barium-oxide- 
ontent glasses, which have a favorable relationship 
between index and dispersion 

The Glass Section of the National Bureau of 
Standards has been engaged in an extensive study of 
Whe refractivity, dispersion, and density in the glass- 
forming region of the barium-borosilicate system [2]. 
lhe phase equilibrium diagram for the system is of 
paramount importance in the interpretation of the 
lata as regards the structure of these glasses. The 
igh silica-boric oxide portion of the system, also, 
ws possible application to the development of im- 
proved enamels and ceramic glazes. Notwithstand- 
ug the fundamental and practical importance of the 
ernary diagram, however, no systematic study of 
he phase equilibrium relationships has been reported. 

The ternary system is of additional interest from 
everal theoretical considerations. It contains one 
inary system (BaO-B,O;) with an immiscibility gap 
ind another (BaO-SiO,) with a portion of the liquidus 
urve approaching that of two-liquid immiscibility. 
u the BaO-SiO, system, BaO - 2SiO, and 2BaO- 35i0, 
how unusual behavior for refractory oxide systems, 
i that they form a complete isomerphous solid 
olution series without a maximum or a minimum. 
As a secondary objective, it was hoped that investi- 
ation of the ternary system might provide addi- 
ional information on these problems. 


s in brackets indicate the literature references at the end of this paper. 
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2. Experimental Procedure 


2.1. Raw Materials 


The starting materials were silica gel (SiO,-nH,O), 
boric acid (H,BO,), and either barium nitrate 
(Ba(NO;),.) or barium carbonate (BaCOQO,). The 
silica gel, 99.99 percent pure on an ignited weight 
basis, was specially prepared by an acid digestion 
process in the Chemistry Division of the National 
Bureau of Standards. The remaining starting ma- 
terials were of reagent quality, meeting ACS specifi- 
cations, 


2.2. Preparation and Analysis of Mixtures 


Because of the volatility of boric acid and the 
extreme hygroscopicity of boric oxide [3, 4], it was 
found advantageous in preparing ternary mixtures 
to add the boric oxide already combined with the 
baria. For this purpose a stock quantity of each of 
the four barium borates, 3BaQO-B,O;, BaO-B,O,, 
BaO-2B,0,;, BaO-4B,0;, was prepared by mixing, 
grinding, sintering or fusing the appropriate amounts 
of boric acid and barium aitrate or carbonate. The 
product was ground to pass a No. 200 mesh screen 
and heated again. The process of grinding and re- 
heating was repeated 2 or 3 times, until a homo- 
geneous product, by petrographic examination, was 
obtained. Barium carbonate was used as the source 
of BaO for preparing 3BaQO-B,O0,; for the other 
borates, barium nitrate was used. The barium 
borates thus prepared were analyzed [4] for BaO and 
B,O;. Previous experience had indicated the ap- 
proximate excess of boric acid required to focmulate 
each compound, and in no case did the actual chemi- 
cal composition vary from the theoretical com- 
pound composition by more than 0.5 percent. No 
attempt was made to adjust compositions to theo- 
retical values, and the analyzed values were applied 
in subsequent calculations 

Stock quantities of the 4 barium. silicates, 
2BaO-Si0,, BaO-SiO,, 2BaO-35iO,, and BaQ-2siQOz, 
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formulated and analyzed com positions 





Tape | lureement between some 
Composition 
‘ | al Formulated * Analyzed 
Bad BO S10 BaO BA 0 
( 2Ra- 310 BaC Oy, SiO 83. 43 16, 37 &3. 53 18. 47 
} RaQ. 2810 BaC Oy, SiO MH. OS 43. 02 55. 04 44. OF 
v4 BaO-2B,05, 2Ba0-38i10 55, 24 +5. 60 + 26 55.11 $5. Sl 9. 38 
»” BaQ4BeO BaQ 2810 17. 25. 92 mw. a8 47.78 25. US 2H. 24 
j BaO-3B,O) 2510 Bal Oy, SiO), Hy BO 58. 30 26. 48 15, 22 S. 48 2H. 24 15, 28 
From compesition and amount 
By analytical determinatior 
Stock supply 
wie batch prepared BeO> added bore at h accounts for a little larger discrepancy, in this case between formulated an ly ze 1 
and of 1 ternary compound, 3BaO-3B,0,-25i0,, P 
were prepared in a similar manner. Barium car- 
bonate exclusively was used as the source of BaO in 
preparing the barium silicates. 
Most of the ternary compositions were prepared 
to lie on joins connecting the compositions of barium 
borates, barium silicates, and the one ternary com- 
pound 3BaQ0-3B,0,;-2Si0,. In making up such 
series, the amounts of conjugate stock compounds 
required to yield 3-g mixtures were calculated, by 
applying the well-known lever principle. The cal- ‘ 
culated amounts, weighed into dry bottles, were 
shaken thoroughly with a mechanical mixer for a si 
“ar . Y 
minimum of | hr, ground with a mechanical grinder ; 
for 1 hr, and then heated for several hours at some l t 
temperature between the solidus and the liquidus. ks 
The products were ground to pass a No. 200 mesh ~ 
sieve and were resintered at a slightly higher temper- “y 
ature, but still below the liquidus. Once again, the 
products were ground and passed through a No. 200 
| TT" ; ; sie Fioure 1 Compositions studied in the systen faQ- Bo 
mesh sieve. us treatment, In most cases, gave a SiO», together with the joins along which most of then 


fairly uniform product, as examined with a polarizing 
microscope. In a few instances, it was necessary to 
heat the material a third time. All heating was done 
in platinum crucibles by using an electrically heated 
furnace. Analysis of representative mixtures pre- 
pared as described showed satisfactory agreement 
between actual and theoretical compositions. Table 
1 gives this agreement for several compounds and 
mixtures. The few compositions not on tie lines 
were made from appropriate quantities of either 2 
already prepared mixtures or 3 stock compounds. 

Analysis of ternary compositions consisted of 
determinations for ignition loss, silica, and baria; 
boric oxide was obtained by difference. Half-gram 
samples were decomposed by 1:1 HCl, evaporated 
almost to dryness, and the moist residue evaporated 
three times with 10-ml portions of methyl alcohol 
to remove the boric oxide. Silica was determined 
by double dehydration in the usual way; baria was 
precipitated and weighed as BaSQ,,. 

Figure 1 shows the 178 compositions investigated, 
together with the joins along which most of them 
were fromulated. It may be noted that many com- 
positions in close proximity to each other were 
formulated on different joins, thus contributing 
statistically to the over-all accuracy of the method. 


formulated 


It. might also be emphasized that the prepared con 
positions always consisted of crystalline or mixtures 
of crystalline and glassy phases, as equilibrium co: 
ditions could be approached more rapidly by melting 
than by erystallization from above the liquidus 


2.3. Apparatus and Method 


‘Phase equilibria relationships in the system we: 
obtained by the well established quenching techniqu 
[5, 6). Charges of about 20 mg were heated in th 
relatively constant temperature zone (hottest regiot 
of a vertical tube, platinum 80 percent-rhodium 2? 
percent wound, resistance furnace [7]. Constan' 
temperature control in this zone to within + 0.3 deg. ( 
was achieved with the use of a newly designe 
self-adjusting, a-c bridge-type controller, describe 
later. 

Temperatures were measured with platinum ver 
sus platinum-rhodium (10%) thermocouples, whic! 


mocouple by the Pyrometry Laboratory of th 
National Bureau of Standards. All temperatur 
values are given on the International Temperatur 
Seale of 1948, based on absolute millivolts [8]. 











were originally calibrated against a standard ther- 
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630.5° and 1,063.0° C (the gold point) the 

ences between the International Temperature 

s of 1948 and 1927 are less than 0.5 deg. C. At 

C. however, the International Seale of 1948 

leg. C less than the scale of 1927; and at 2,000° C 

ference is 6.4 deg. C [9]. In reviewing the litera- 

ff high-temperature phase studies, it is found 

many investigators go into great detail 

the accuracy of thermocouple calibrations, 

ever specify the temperature scale used. It 

it be safely assumed that all temperature values 

January 1, 1948, are reported on the Inter- 

mal Temperature Seale of 1948, even though on 

day, by international agreement, the change 

f international to absolute units became effec- 

The changeover at the National Bureau of 

Siandards, for example, did not become effective 

January 1, 1949. Since 1914 the Geophysical 

Laboratory of the Carnegie Institution of Washing- 

nas used a temperature scale based on the work 

Day and Sosman. The scale is in good agreement 

vith the International Temperature Scale of 1948, 
oO about 1.550° & 

Che thermocouple calibrations were checked peri- 


con- 


odically against the melting points of barium 
tetraborate (889° C) [4], the ternary compound, 
\BaO-3 BlOs-2810, (1,009° C), gold (1,063°C), barium 
borate (1,105° C) [4] and barium disilicate (1,416 
( 


Quenched samples were crushed and examined 
vith the polarizing microscope, using the standard 
nersion media methods. Refractive indices meas- 
this method are accurate to within +0.003 


wre dl by 


X-ray diffraction powder patterns (with copper K 
alpha radiation) were made for a number of com- 
and mixtures. A commercial type Geiger 
counter X-ray spectrometer equipped with a syn- 
chronous motor-driven scanning unit and attached 
to an electronic high-speed recorder was used for 


po inds 


this purpose 


2.4 Adjusting, Bridge-type Controller 
\ self-adjusting, bridge-type controller suitable 
for precise temperature control of a quenching furnace 
and assembled from 


was designed commercially 
available components by F. A. Mauer of the Na- 
tional Bureau of Standards. Like the electric- 


furnace thermostat developed by H. S. Roberts at 
the Geophysical Laboratory in 1925 [7], the controller 
operates on the change of resistance that accompa- 
nies a change in temperature of the furnace winding. 
The contacting galvanometer has been replaced by 
a servo amplifier (A, fig. 2), which provides power 
to operate a reversible, two-phase motor, M. As 
the amplifier operates on a 60-cycle a-c signal, the 
power for the furnace is adjusted by means of a 
variable autotransformer. This transformer is posi- 
tioned by the reversible motor and provides the 
power required to maintain the furnace winding at 
the temperature prescribed by the resistance setting. 
lhe temperature is maintained constant, of course, 





Ir 


accordance with the International Temperature Scale of 1048. 
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Fiaurre 2 Schemat diagram of self-adjusting, bridge-ty pe 


controle for high tem perature Jurnace 





with the 


nq ¢ 


quen h ng 


controller together 


FiaurE 3 
rnace and the temperature measu 


Te m perature 


aipment 


only if volatilization of the winding is negligible 

Figure 2 shows a schematic diagram of the con- 
troller, and figure 3 shows the controller, togethe: 
with the quenching furnace and the temperature 
measuring equipment. 

The controller offers a combination of features 
that have not been available in any one instrument 
These principal advantages are as follows 

1. Control of temperature to within +0.3 deg. C 

2. Elimination of process lags by use of the furnace 
winding as the sensing element. 

3. Use of alternating current to heat the furnace 
This feature is important because a-c power is 


readily available, permits the use of a variable 
adjust 


autotransformer to the voltage without 











wasting power, and also permits amplification of the 
signal from the bridge 

4. Incorporation of proportional speed floating 
control with automatic reset. When the temperature 
departs only slightly from the control point, the 
autotransformer moves at a rate proportional to the 
departure. In the case of a large unbalance, the 
motor turns at an essentially constant speed until 
balance is reached. Because the motor rotates as 
long as the temperature remains different from the 
required temperature, the variac automatically 
assumes a new position for each temperature setting 
and even compensates for change in the line voltage 
of up to 20 percent. As long as power requirements 
and line voltage remain constant, oscillations of the 
balancing motor are barely perceptible. Conse- 
quently, even if the controller is switched off, the 
temperature will remain within a few degrees of the 
control point. This characteristic simplifies the 
problem of controlling the furnace in a desired 
temperature range so that it can be ready for use on 
the following day. 

5. Simplification of operating procedure. A single 
knob is used for resetting the control point. Because 
of the reset action described above, the controller 
will position the variable autotransformer without 
further attention, to provide the power requirement 
for a new setting of the balance point. The control 
point setting can be varied continuously by means of 
a 10-turn precision potentiometer. 

It should be noted that since an a-c bridge circuit 
is used, an impedance balance must be effected if 
maximum precision is desired. The bridge output 
can be recorded on an oscilloscope and capacitors 
added to the high-impedance arms of the bridge 
until a true null can be obtained [10]. 

The principal maintenance problem is presented 
by the variable autotransformer brushes. During 
operation, the “hunting” of the controller keeps the 
brushes moving back and forth between adjacent 
windings. Experience indicates that after about 6 
months of operation, arcing may occur between the 
brushes and the slightly roughened windings. In 
this case, the brushes are cleaned or replaced, and the 
windings are cleaned according to the manufacturer’s 
instructions. 


3. Limiting Binary Systems 
3.1. System BaO-B,O, 


The binary system BaO-B,O, was reported by 
Levin and McCurdie [4]. It is distinguished by a 
region of two-liquid immiscibility extending from 
about 70 percent of B,O, (30 percent of BaO) to 
almost pure B,O;. Extrapolation from ternary data 
to the binary indicated that the reported binary 
liquidus temperatures were uniformly low by about 
10° C. New determinations on the binary com- 
pounds confirmed this conclusion. A corrected dia- 
gram is reproduced as figure 4. 


3.2. System B,O,-SiO, 


Only limited information is available about the 
system B,O,-SiO,, and no liquidus values are known, 
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Liquidus values not changed to International temperature scale of 1948 


as no mixture in the system has been crystallized 
Morey [11] was unable to crystallize any mixture i 
the system, either from a dry mix or by hydrotherma 
treatment. Greig [12], and Cousin and Turner [13 
believe that the oxides as liquids are completely 
miscible in all proportions. The latter authors, o1 
the basis of thermal expansion and densit y-composi- 
tion curves, could find no indication of compound 
formation. N. K. Dertev and T. H. Khudyakova 
[14], however, believe that several compounds (for 
mulas unspecified) are possible. 


3.3. System BaO-SiO, 


The binary system BaO-SiQ,, as shown in figure 5 
was investigated originally by Eskola [15] and mod- 
ified later by Greig [12]. It shows two regions o! 
particular interest. In the first, dibarium trisilicat: 
(2Ba0.-3Si0,) and barium disilicate (BaO-25i0, 
form a complete solid solution series. The liquidus 
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«der pattern of crystals grown at 1,100° C showed very slight irregular 
1cings from pure BaQ-2Si0:. No detectable change in optical 
ition with BaO-28i0, and with barium borates, confirmed by optical 
st 
vattern of primery phese crystals showed a maximum shift of 
gs from pure 2BaO-3S8i0> No detectable change in optical 





juench growths and barium platinates 


ies of the BaO-SiQ, system. Melting points, 
tical properties, and X-ray diffraction data for all 
ompounds in the system were determined; and a 
microscopic and X-ray study of the BaO-2si0, 
2BaO-3Si0, solid solution area was made. 

lhe melting points and optical properties of the 
compounds were found to be, in general, in close 
agreement with the reported values [15, 16]. Eskola 
15] reported the melting point of dibarium silicate 
2BaO0-SiO,) as being above that of platinum. By 
a method of optical pyrometry [19j, it was possible 
to estimate the melting point as 1,815 +25° C. 
Microscopic examination of 2BaQO-SiO, that had 
been heated above 1,800° C gave some new optical 
data. The material appeared as irregular grains 
showing good cleavage parallel to the plane contain- 
ing the alpha and beta vibration directions; it was 
biaxial positive with a very small 2V (5 to 10°); 

1.800, B=1.80+, y= 1.826. 

Table 2 gives a comparison of the X-ray powder 
pattern diffraction data obtained in this study with 
that of Austin’s [20]. Although agreement between 
the two sets of data is reasonably good, enough 
differences exist to warrant recording. It can be 
from table 2 that minor differences exist in 
interplanar spacings and relative intensities between 
corresponding patterns. The present study reveals 
additional interplanar spacings for most of the com- 
pounds, in particular for BaO-2SiO,. Perhaps the 
most notable difference in the corresponding pairs 
of patterns is the lack of agreement of the three 
major peaks (identified as (1), (2), and (3) in table 
2). Only the numbers (1) and (2) peaks for 
2BaO-SiO, and the number (1) peaks for BaO-SiO, 
and 2BaO-3SiO, agree between the two sets of data. 
(he remainder of the peaks are either interchanged 

order of intensity or are entirely different. The 
discrepancies in order of intensity of the three major 
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in the system BaQO- B,O,-SiO,—¢ tinued 
2510) | na f nu 
Phase ‘ 
n of 
quenched 
glass 
256° C Primary Cot he (ilas 
1. (41 All 
Moderate w 
4s All 
Rare x 
1. (47 All 
Rare x 
1. O48 All 
Rare . 
1. 648 All 
Rare x 
Moderate Ba0. B,O 


¢ Liquidus greater than this temperature 

‘ Temperature of a ternary invariant point 

« All glass devitrified to form quench growths 

+ Glass and quench growths 

lo check liquidus of binary immiscibility re 

i Most difficult of the primary phase regions 
sccurate; tridymite only form identified 

t Composition of 3Ba0-3B,0 


gen 
studied and, consequently, the le 


2810», 


peaks for corresponding sets of data are due probably 
to orientation. In this study, an attempt was made 
to minimize orientation in the sample during speci- 
men mounting. The finely ground material (passing 
a No. 200 mesh screen) was neither packed nor 
scraped at any stage of preparation but was sprinkled 
freely over a microscope slide (which served as a 
specimen holder) that had been smeared with a thin 
laver of petrolatum. 

Study of the BaO-SiO,-2BaO-3SiO, solid solution 
area was based on the following four compositions 
100 percent of BaO-2SiO, (0 percent of 2BaO-3Si0,), 
70.8 percent of BaO-2SiO, (29.2 percent of 2BaO. 
35i0,), 37.5 percent of BaO-2SiO, (62.5 percent of 
2BaO0-35i0,), and 0 percent of BaO-2SiO, (100 
percent of 2BaQO-.35i02 Results of the optical 
examination agreed with Eskola’s study [15;. The 
X-ray data, however, did not conclusively show a 
solid solution series. The powder pattern for the 
composition containing 70.8 percent of BaO-SiO, 
(29.2 percent of 2BaO-3SiO,) showed several double 
peaks which indicated the presence of two phases 
Even after an additional melting and grinding of the 
mixture, these double peaks persisted. It is believed 
that single crystal work, which would permit the 
determination of unit cell sizes, would be the most 
suitable method of studyiag the solid solution area 
by means of X-rays. 


4. Experimental Results and Discussion 
4.1. General Statement 
Table 3 gives the essential quenching data for the 


178 compositions studied, grouped according to 
primary fields of crystallization. ‘The compositions 


within each primary field are arranged according to 
The locations of the 


increasing percentage of BaO. 
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compositions studied are shown in figure 1. The 
data from table 3 are plotted in figures 6 and 7. 
Figure 6 shows the ternary diagram for the system 
with the temperatures of the invariant points listed 
alongside. The heavy lines, solid and dashed, rep- 
resent the phase boundaries and delineate the pri- 
mary phase areas. The arrows indicate the direction 
of falling temperature, in the conventional mannet 
Isotherms are drawn generally at 100-deg intervals 
but in some cases at 50- or 25-deg intervals. The 
temperatures for the polymorphic inversions of silica 
such as quartz to tridymite and tridymite to cristo- 
balite, could not be obtained experimentally, but for 
the case of the tridymite-cristobalite inversion the 
assumption was made that the inversion temperature 
is the same as for pure silica, namely, 1,470° C 
Little data could be obtained in the regions very rich 
in BaO, for example, in the 2BaO-SiO, and BaO 
primary phase areas, (1) because of the reaction of the 
samples with platinum at high temperatures to form 
barium platinates, (2) because of the extremely rapid 
rate of devitrification of samples during quenching 
and (3) because of the rapid hydration of samples 
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Experimental difficulties due 

_ it is believed, to the high viscosity of the 
were encountered, likewise, in studying most 
two-liquid immiscibility region and the silica 
vy phase area. It should be noted that such 
s were encountered as very fine crystals difficult 
ntify, evidences of nonequilibrium coaditions, 
e to obtain clear-cut separation of the immiscible 
is, and inability over most of the silica field to 
mine liquidus temperatures. This portion of 
agram, therefore, is of a tentative nature, but it 
sught that the relationships as shown are essen- 
correct. 

ch of the eight binary compounds shows a 

ary phase region in the ternary diagram (fig. 6). 

3Ba0-3B,0,-25i0,, is 


que nec ‘hing. 


‘) one ternary compound, 

found. The system contains 5 ternary eutectics, 
lettered A, B, C, D, E; five ternary invariant points 
1 itectics (reaction points lettered F, G, H, I, J 
and 2 quintuple points, lettered K and L, at which 2 
liquid phases are in equilibrium with 2 crystalline 


phases and vapor. With the exception of point J, 
temperatures of these invariant points range between 
s and 980° C. The invariant point J is under 
1.370° C, the value of the binary eutectic between 
BaO and 3BaO0-B,0,. However, as point J is located 
near the binary eutectic, it is to be expected that its 
liquidus value is not much lower. One other in- 
variant point, M, representing the equilibrium be- 
tween silica, boric oxide, and BaO-4B,0, is postulated 
(21) to lie somewhere in the narrow region between 
the B,O,;-SiQ, binary and the boundary between the 
two immiscible liquids and one-liquid regions. The 
most peominent feature of the diagram is the large 
extent of the area of two immiscible liquids, which 
covers approximately 50 percent of the diagram. It 
is also interesting to note that the BaO-4B,0, and 
SiO, primary phase areas include almost 70 percent 
of the diagram. 

Figure 7 shows the eleven composition 
formed in the svstem by the Alkemade lines. For 
the purpose of clarity, isotherms are not shown. 
Three secondary binary systems are evident: 
BaO-SiO,-BaO-B,O0;,, 3BaO0-3B.0,-2Si0,.-BaO-B,Os, 
and 3BaO-3B,03-2Si0,-BaO-2B,0; (see also fig. 8, 
D, E, F). It cannot be stated with assurance 
whether or not the system 3BaQ-3B,0;-2SiO,- 
BaQO-2Si0, is a binary because the exact nature and 
extent of the solid solution field is not known (see fig. 
8, A). The system 3Ba0-3B,0 3-25i0.-BaQO-B.O;- 
BaO-2B,0;, is seen to be a secondary ternary system. 
Crystallization paths for compositions in that por- 
tion of the baria-borosilicate system where liquid 
immiscibility does not occur are, in general, straight- 
forward and familiar types. Typical crystallization 
+ hs in a two-liquid area are discussed by Flint and 


ells [3]. 
4.2. Stability Fields of the Binary Compounds 
a. BaO.2SiO, 


triangles 


Barium disilicate forms solid solutions with 


2BaQ-35i0, and probably with the barium borates, 








the partial limits of these solid solutioas are shown 
in figure 6. Compositions studied in the BaO-2SiO, 


primary field are given in table 3, part A. Liquidus 
temperatures vary from the melting point of 
|} BaO-2SiO, (1,418° C)* to a minimum of 950° C 


(point E). Optical properties of the primary phase 
crystals within the limits of experimental measure- 
ments did not differ from those of pure BaO-2SiQ,. 
X-ray diffraction powder pattern data of primary 
phase crystals grown at 1,100° C for two composi- 
tions (11 and 13 of table 3) showed a very slight and 
irregular shift of d spacings from the pure BaO-2SiQy. 
The shift averaged less than 0.01 A, which is about 
the expected limit of accuracy of the experimental 
conditions. When liquidus temperatures along the 
join BaO-2Si0,-3BaQ0-3B,0;-25iO, are plotted (fig. 
8, A), a distinct and experimentally significant 
discontinuity is found in the liquidus curve at about 
1,250° C. At this break, the minimum deviation of 
the curve from one which would show no breaks 
(that is, a straight line) is 50° C, a value undoubtedly 
above the limit of experimental error. It must be 
concluded, therefore, that the break in the curve is 
is indicative of a relatively sharp phase change. It 
would seem likely that the break represents the 
transition temperature vetween two closely related 
solid solutio 1 phases, but the nature of these phases 
is not known. 


b. 2BaO- 3SiO, 
Dibarium trisilicate forms solid solutions with 
BaO-2Si0, and with the barium borates. The par- 


tial limits of the solid solutions are shown in figure 6. 
Compositions studied in the 2BaQ-3SiO, primary 
field are given in table 3, B. 

Primary phase crystals from compositions 21, 22, 
29, 30, 33, 37 (table 3, B) gave distinct evidence of 
being solid solution crystals, both by optical and 
X-ray methods. No detectable differences in optical 
properties from pure 2BaQ-3SiO, were found for the 
primary phase crystals obtained from the remaining 
compositions. However, cry stals from the composi- 
tions 25, 26, 38, and 42 showed small and irregular 
shifts in d spacings, to a maximum extent of 0.02 A. 
X-ray data were not obtained for all compositions in 
the field, but it seems from the limited data available 
that all or most of the compositions would show 
similar slight differences. 

When the liquidus temperatures are plotted for 


the join 3BaO0-3B,0,-2Si0,-2BaO-3Si0,, a discon- 
tinuity is found in the liquidus curve at about 
1,250° C (fig. 8, B), or the same temperature as the 


join. The 
BaO-B,O,-2BaO-35i0, 
2Ba0-3Si0, 


break in the 3BaQ0-3B,0,-2Si0,-BaO-2Si0, 
liquidus curve for the join 


(fig. 8, C) similarly shows a break in the 


| solid solution field, but it is not as pronounced and 


occurs at a lower temperature (approximately 1,125 
C). The latter join is more nearly parallel to the 
isotherms than the other two joins mentioned, and, 
consequently, the location of the break may be less 
accurate. 


+ Temperature values on the BaO -SiO) binary have not been converted to 
the International Temperature Scale of 1948 
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study of the system BaQO-Al,O;-SiQ, 
s [16] found varying indices of refraction for 
\positions within the field surrounding the 
sition 2BaO-3SiO,; and, therefore, he proposed 
BaO-3SiO, is not a true compound but merely 
ilar composition at or near the maximum in a 
solution series. X-ray powder pattern data 
ed to confirm his conclusion, although the 
ons in d values were not as pronounced as the 
ons in indices of refraction. In the present 
the reverse situation is true, namely, that 
measurements In Many instances appear to 
very slight differences that are not revealed by 
methods. The present study, likewise, does 
melusively prove the existence of a 2BaO-3Si€ ), 
iry phase field of unvarying composition and 


1s 


rties 
question as to whether or not dibarium tri- 
te is a compound becomes especially important 
drawing the composition triangle (fig. 7) 
vusly, if dibarium trisilicate is not a compound, 
-emade line can originate from its composition 
Finally, the method of eutectic matrix com- 


_ 


ven pertinent ternary compositions, shown by 
he dots in figure 7, were selected for study. The 

ilated compositions were melted and _ stirred 
wve the liquidus in order to obtain homogeneity. 
he melts were then crystallized by subjecting them 
slow cooling schedule. The phases in the de- 
fied samples were identified by petrographic 
examination. 

\lost of the samples contained varying amounts of 

ench growths and small amounts of glass, indicat- 
ing that equilibrium conditions and complete de- 

itrification were not achieved. In most cases only 
two crystalline phases were identified. This may be 
explained on the basis that devitrification of the 
solid solution erystals consumed all the available 
iquid before a third phase crystallized. In the case 
of one composition (a, table 4) for which no solid- 
solution phase appeared, a different explanation is 
necessary, as, for example, too small a percentage of 
fine-grained material intermixed with the quench 
growths of another phase. 

Figure 9 and table 4 show the location of the com- 
positions studied, together with the phases identified. 
Compositions 6, c, and e after devitrification showed 
three phases: 3BaO0-3B,0,-2Si0,, BaO-B.O,, and a 





Phases identified in seven 


whether 


Composi ( eee weight Phase 
Desig- tion num percen 
nation ber in 

table 3 BaO BO SIO, {BaO-3Be 





devitrified glasses, in ‘an attempt to determine by the method of eutectic matris 
Baw) BSiO, 18% compound or merely a party lar com position 


solid solution with optical properties closely resem- 
bling 2BaO-3Si0, Only 3BaO0-3B,0,-2810, and 
BaO-B,O, were found in composition a. The data 
from the four mentioned compositions are consistent 


with the compatability triangles as drawn. If the 
correct Alkemade line’ were really the — join 
3Ba0-3B,0,-2S810,-BaO-SiO,, u possibility if 


2BaO0-3Si0, is only a composition in a partial solid 
solution series between BaO.2SiO, and BaO.-SiO,, a 
BaO-B.O, phase should not have been obtained for 
compositions 6 and ¢, but rather some indication of 
BaO-SiO, Furthermore, composition g gave no 
indication of a 3BaO-3B,0,-2SiO, phase whereas a, 
only a few percent away, showed a moderate amount 
of 3Ba0-3B,0,-2Si0, It seems reasonable to sup- 
pose that an Alkemade line, as drawn, passes be- 
tween the two compositions. Compositions 6 and 
e both showed BaO-B,O,, whereas compositions d 


and f, nearby, did not. The only join that satis- 
fies this result connects 3BaO-3B,0,-2Si0, and 
2BaO0-3Si0, Because the composition 2BaQ-3Si0, 


seems to represent the only logical point for the 
origination of two Alkemade lines, it must be con- 
cluded that 2BaO-3SiO, is a compound, regardless 
of the extent of the solid-solution region 
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c. BaO.SiO, 
Compositions determining the barium silicate 
primary phase area (fig. 6) are listed in table 3, 


part C. Optical properties of the BaO-SiO, crystals 
agreed well with the recorded values [15, 16, 20). 
The lowest temperature for any liquidus in the 
field is 875 +5° C, at the ternary eutectic C. Barium 
silicate and barium borate form a true binary sys- 
tem, with a eutectic at about 940° C and 45 percent, 
by weight, BaO-SiO, (fig. 8, D 


d. 2BaO.SiO, and BaO 


Information on the field of dibarium silicate could 


be obtained from only seven compositions (table 3, 


part D). Six of these compositions were located 
along the 2BaQO-Si0,-BaO-SiO, boundary. The 
seventh (No. 70), on about the 79 percent BaO 


isopleth, formed quench growths and barium plati- 
nates to the extent that the liquidus value could 
not be determined; and such was the case for all 
compositions in the field richer in BaO Dashed 
isotherms were drawn by extrapolation from known 
boundary curves. 

For similar reasons, together with the high tem- 
peratures involved, no liquidus values could be 
obtained in the BaO primary phase region. Because 
of the reaction of BaO at elevated temperatures with 
the platinum of a thermocouple, differential thermal 


analysis cannot be applied, as was discovered in the | 


study of the BaO-B,O, binary system [4]. It is ap- 
parent, however, that the BaO primary field must be 
limited to a small area in the corner of the diagram. 


e. 3Ba0.B,O 


Liquidus determinations in the tribarium borate 
field presented most of the same difficulties as in the 
2BaO.-SiO, and BaO areas; however, the three com- 
positions studied (table 3, part E) served to delineate 
the tribarium borate field, which is the smallest in 
the system. Composition 71 (liquidus 908° C) was 
close to the ternary eutectic C (fig. 6), and a sample 
quenched from 875° C showed three crystalline 
phases (3Ba0-B,O0,, BaO-B,O,, and BaO-SiQ,) in 
equilibrium with small amount of glass. A ternary 


reaction point is postulated at J (figs. 6 and 7), where | 


the crystalline phases 3BaO-B,0,, BaO and 
2Ba0.-2i0, would be in equilibrium with glass. The 
invariant point J might be expected to be near in 
both temperature and composition to the binary 
eutectic between 3BaO0-B,O, and BaO at 1,370° C. 


f. BaO.B,O, 


Eighteen compositions (table 3, F) were studied 
in the primary phase field of barium borate, and no 
experimental difficulties or unusual features were 
encountered. Four ternary invariant points C, D, 
H, and B (figs. 6 and 7) ranging in temperature from 
875° to 925° C were located on its boundaries. 
Barium borate forms true binary systems with 


both BaO-SiO, and 3Ba0-3B,0,-2Si0, (fig. 8, D 
and E), 











g. BaO.2B.0 


The 12 compositions studied in the barium dil 
field are given in table 3, G. The compound { spp, 
a true binary system with 3BaO0-3B,0,-2Si0O. ({, 
8, F) and accordingly a temperature maximy 9 | 
shown on the boundary curve BF of figure 6. 


h. BaO.4B,0; 


Compositions for this field are listed in table 3. } 
Barium tetraborate occupies a primary phase revio; 
second in extent only to that of silica. It is interesy. 
ing to note, also, that less than 5 percent of the field | 
lies in a one-liquid region, whereas more than 9: 
percent is in the two-liquid immiscible area (fig. § 
As the percentage of silica increased, determinations 
of liquidus temperatures became successively mor 
difficult because of sluggish crystallization of th 
viscous melts. The boundary between the two. 
liquid and the one-liquid region adjacent to thy 
B,O,-SiO, binary was not determined and is show 
as a dashed line (figs. 6 and 7). That this bounda; 
must lie within several percent of the B,O,-SiO 
binary is deduced from the limited data on the in- 
dices of refraction of the low index conjugate glass 
(see ‘‘n of quenched glass’’, table 3, H). The values 
are approximately the same or less than those fo; 
pure SiO, or B,O;, both of which are about 1.46 [24 


i. B,O, 


If the compatability triangle formed by BaO-4B,0 
BOs, and SiO, (fig. 7) is correct, no additional com- 
pounds showing a stable primary phase area can ly 
formed within its boundaries. It is then necessary 
to postulate an invariant point M located in th 
narrow one-liquid region adjacent to the B,O,-Si0, 
binary, where the three primary phases coexist ir 
equilibrium. In this event boric oxide would 
possess a very narrow primary field extending from 
the invariant point to pure boric oxide. Because of 
experimental difficulties, no attempt was made to 
substantiate this conclusion 


j. SiO 


The silica primary phase area (fig. 6) occupies 
approximately 40 percent of the phase diagram, and 
the silica liquidus surface is distorted by the occur- 
rence of liquid immiscibility. The field consists of 
a large portion in the two-liquid area limited o 
three sides by a portion in the one-liquid area. A 
narrow strip of one-liquid area is present along the 
B,O,-SiO, binary and extends to the postulated 
eutectic at M. Another portion of the one-liquid 
area extends along the BaO-SiO, binary to the 
eutectic at 1,374° C between BaO.2SiO, and SiO,, 
and then penetrates the interior of the diagram 
adjacent to the primary phase fields of BaO-2Si0,, 
3BaO0-3B,0,-2Si0,, and BaO.2B,0,. The invariant 
point A representing the equilibrium between SiO, 
BaO-2B.0, and BaO-4B.0, is at 810° C. That the 
silica field must extend to the BaO-4B,0, field was 
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vident from the fact that SiO, and BaO-.4B.O 
he only primary phases identified in the two- 
area, and fields entering a two-liquid area 


persist until liquid miscibility is again en- 
red 
tematic investigation of the silica field was 


ered by the extremely high viscosities of all the 
. and the attendant slow rate of crystal growth 
pproach to equilibrium. Actual liquidus values 
not be obtained for a number of compositions 
117, 119, 122, 123, 128 table 3, I); but phase 
fication in these cases was useful in approxi- 
¢ the limits of liquid immiscibility. In areas 
cent to the primary phase boundaries, it was not 
hy possible to determine the exact temperature 
of appearance or disappearance of the silica phase. 
Consequently, the silica area of the diagram is less 
irate than the rest both as regards to temperature 

and boundaries 
idvymite was the only form of silica identified. 
Even at temperatures as low as 823° C (No. 125 of 
3, 1), or about 44 deg below the inversion 
perature of pure tridymite to quartz, no quartz 
formed in any of the mixtures. In the study of 
the system K,O-MgO-SiO, [25], Roedder also found 
tridymite persisting metastably, to a eutectic point 
at 715 +10° C. He concluded that a quartz 
liquidus stable at a higher temperature presumably 
existed, but because of the vi nature of the 
iss, excessive time would be required to demon- 


viscous 


strate the point 


4.3. Stability Field of the Ternary Compound 
3BaO-3B,0,-2SiO, 


This was the only ternary compound found in the 
system, and 37 compositions within its primary 
field were studied (table 3, J). The field occupies a 
large portion of the glass-forming region of the 
system, and it was observed qualitatively that rates 
of devitrification (or of erystal growth) with the 
exception of compositions in the silica field were 
much slower than for compositions in the other 
primary phase areas 

The compound 3Ba0-3B,03-2SiO, melts congru- 
ently at 1,009° C. As seen from the spacing of the 
isotherms (fig. 6), its primary field is extremely 
flat, indicative of a high degree of dissociation of the 
compound in the liquid state [26]. The compound 
crystallizes as irregular, elongated, and _ fiber-like 
grains with positive sign of elongation. It is biaxial 
negative: a=1.669, B=1.696, y=1.702 C), 
with a moderate optic axial angle of about 35°. 
The refractive index of the glass is 1.633. When the 
compound is formed by 


(25 


synthetic twinning is prominent. Data from the 
\-rav diffraction powder patterns are given = in 
table 5. 


Five ternary invariant points (B, H, G, E, F, fig 
6) are found along the boundaries of the 3BaQO. 
1B.04-28i0, field. The lowest of these invariant 
points is at 825°+10° C (point F). The ternary 
compound forms true binary systems with BaO-B,O, 


devitrifving glass, poly- | 


53 





TABLE 5 Y-ray di ffract on data fo the om pound 
3BaO0-3B,05-28i0O2, showing inte planar spacings (d) and 
their relative intensities (] 
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and with BaQO-2B,O, (fig. 8, E, F). As discussed 
earlier, it is not believed that 3BaO-3B,0,-2Si0. 
forms a true binary system with BaO-2SiO,, because 
the nature of the solid solution area is not known 
(fig. 8, A). 

4.4. Two-Liquid Area 


The approximate shape of the two-liquid area at 
the liquidus is shown in figure 6. An interesting 
feature of the immiscibility area is the marked 
increase in concentration of BaO in the less siliceous 
liquids, proceeding toward the interior of the dia- 
gram from the baria-boric oxide binary. ‘The maxi- 
mum content of BaO in the ternary liquids is about 
43 percent, which is an increase of 11 percent over 
the BaO content of the corresponding binary liquid 
A similar observation was made in the lime-boric 
oxide-silica system [3, 21] and in the ZnO-B,O,-SiO, 
system [27]; but in these systems both the CaQO-SiQ, 
and the CaO-B,O, binaries and the ZnO-SiO, and 
ZnO-B,O, binaries, respectively, exhibit large immis 
cible regions. In the present system only the 
BaO-B,O;, binary possesses a liquid immiscibilits 
region, but the flat liquidus portion of the BaO-SiO, 
system indicating an approach to two-liquid immis- 
cibility as discussed earlier (fig. 5) is probably a 
factor contributing to the wide extent of liquid 
immiscibility in the ternary system. In a study of 
the effect of lanthanum and barium in glass-forming 
systems, Brewster and his coworkers [1] present 
some very limited data on immiscibility in the 
system BaO-B,O,-La,O,. Their data indicate that 
in this system, the addition of La,OQ, to the binary 
system BaO-B,O, causes a decrease in extent of 
liquid immiscibility 

Geller and Bunting [28] in an investigation of the 
system PbO-B,O;-SiO, found a large region of two- 
liquid immiscibility similar to that in the CaO-B,O,- 
SiO, system. The PbO-B,O,; binary has an immis- 
cibility gap [29], but none has been reported in the 
PbO-SiO, binary [30, 31]. The immiscibility area, 
nevertheless, shows an increase in concentration of 

















PbO in the less siliceous liquids proceeding toward 
the interior of the diagram from the PbO-B,O, binary 
The PbO concentration varies from about 45 percent 
on the PbO-B,O, binary, to a maximum of about 
60 percent. Inspection of the PbO-SiO, diagram 
showed that the liquidus curve had not been deter- 
mined beyond about 35 percent, by weight, of SiQy. 
It is postulated, therefore, from the shape of the 
ternary immiscibility boundary that the liquidus 
curve mm the unexplored area of the PbO-SiO, binary 
might be S-shaped as in the BaO-SiQ, system. 

A study was made of the compositions of conju- 
gate liquids formed within the immiscible area by 
melting of two compositions, one in the BaO-4B,0, 
and the other in the SiO, primary field. The results 
are listed in table 6 and shown in a portion of figure 
7. In figure 7, the solid dots in the BaO-4B,0, and 
SiO, fields represent the two initial compositions; 
and the open circles connected by straight lines, the 
analyzed values of the separated compositions. 

The initial compositions were prepared from the 
end members, barium carbonate, silica gel and boric 
acid, by the described method of repeated grinding 
and partial fusion. These two compositions were 
not analyzed at the start, but the two separated 
compositions, regardless of the degree of separation, 
and the associated initial one must lie on a straight 
lind provided there has been no change on the initial 
composition during the heat treatment. Figure 7 
shows that the three associated compositions do lie 
essentially on a straight line; and it is reasonable to 
assume that composition number 118 (in the silica 
field) was altered slightly during its preparation 
In no case was separation of the conjugate immiscible 
liquids complete, but the intersection of the line 
passing through them with the boundaries of the 
immiscible area is a tie line for the particular tem- 
perature of separation. 

The quintuple line KL is the boundary separating 
the BaO-4B,0, and SiO, primary fields. Any point 
on the quintuple line represents the equilibrium be- 
tween two liquids of compositions K and L, two solids 
BaO-4B,0, and SiO,, and vapor 

As separation of the immiscible liquids could not 
be achieved by gravity at the liquidus temperature 
or slightly above, the actual separations were made 
well above the liquidus. For composition number 
104 (fig. 7 and table 6) separation was tried at two 
temperatures, 925° and 1,250° C, in an attempt to 
obtain information on the effect of temperature on 


the extent of two-liquid immiscibility. Because o| 
incomplete separation, it was not possible to dray 
any conclusions from the lengths of the two tie lines 
Inspection of the data, however, in two othe: 
spects was revealing: (1) The tie lines for the tw 
temperatures of separation for composition 104 

7) showed a small but significant angle between them 
This rotation of the tie line with increasing temper 
ature is positive indication that the compositions of 
the conjugate liquids are changing with temperatur 
[32]. (2) The indices of refraction, measured with 
the polarizing microscope, of clear fragments of thy 
low silica (high index) conjugate glasses were as 
follows: 1.542 +0.002 and 1.540 +0.002 (25° C 
for the glasses separated as liquids at 925° and 1,250 
C, respectively. Extrapolations from figure 10 show- 
ing the indices of refraction of the quenched glasses 
indicate that a decrease in index of refraction of 
0.002 corresponds to a decrease in BaO content of 
about 2 percent. Measurements of the high silica 
conjugate glasses were more difficult and less exact 
but the index of refraction of the glass separated at 
1.250° C was definitely higher (between 1.45 and 
1.46) than the glass separated at 925° C (just under 
1.45 It seems, therefore, from these data that th: 
two-liquid area decreases in extent with increasing 
temperature. 

A few preliminary experiments were performed 
using the high-temperature centrifuge developed by 
T. F. Newkirk and Fred Ordway * [33]. Much better 
separations than by gravity were obtained, and this 
approach offers new promise in the study of liquid 
immiscibility. Their equipment, however, is suit- 
able only for small samples, and micromethods of 
analysis would have to be applied. 

Recent studies on several selected compositions |o- 
cated near the boundary of the BaO rich immiscible 
liquid have shown that it is possible to quench these 
compositions as two glasses at one temperature and 
as one glass at some higher temperature, presumably 
where complete miscibility is encountered. Work 
is in progress with the view of determining the shape 
of the immiscibility dome. 


4.5 Insofract diagram for quenched glasses 


shows an isofract diagram for the 
It can be seen that barium oxide 


Figure 10 
quenched glasses. 





* Associated with the Portland Cement Association Fellowship at the Nationa 
Bureau of Standards. 
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study, in which no attempt was made to relieve the 
internal stresses and strains caused by quenching 
; 5. Summary 
‘ 
The equilibrium diagram for the system BaO- 
B.O;-SiO, was constructed from data, obtained essen- 
f tially by the quenching method, on 178 ternary com- 
{ positions. All the binary compounds showed pri- 
: mary fields in the ternary system. One new com- 
eo pound, 3BaO-3B,05-2S5i0,, melting at 1,009° C and 
possessing an extremely flat primary field was found, 
shies . and its optical and X-ray properties described. Table 
sc et | “iy 7 lists the compositions and liquidus temperatures 
ee ‘ of the ternary invariant points together with the 
enee : .” crystalline phases in equilibrium with liquid 
> | ‘ , The BaO-4B,0; and SiO, primary phase areas 
wet ' q include almost 70 percent of the diagram, and ove! 
/ ‘ these areas exists a large region of two-liquid immisci- 
hom 8 bility. On the basis of rotation of a tie line in the 
per 154 immiscibility region at two temperatures and of a 
is of é “3.02 Change in index of refraction of the conjugate glasses 
on we “# at the two temperatures, it was concluded that the im- 
with fa 10. Indices of refraction of quenched glasses in the | ™iscibility gap decreases with increasing temperature 
the | system BaO-—B,O;-SiO Applying the method of eutectic matrix compound 
. @ ; composition showed that 2BaQO-3SiO, is probably a 
( its the greatest change in refractive index and | compound and not merely a particular composition, 
tT hat the increase in refractive index for a given addi- | even though the exact nature and extent of the 
aie of BaO becomes greater as the BaO corner is | 2BaO-3SiO, solid-solution field could not be ascer- 
isses [an oached. Of particular interest is the convex | tained. The solid-solution region between 2BaO- 
n of ature of the isofracts, with the BaO vertex asa com- | 3510, and Baf )- 225i ), extends into the interior of the 
t of } mon center of curvature. Thus, indices of glasses | diagram and was found to be complex 
ilica : on a BaO isopleth will show a maximum at An isofract diagram for the quenched glasses shows 
act P about the center of the isopleth. These findings | that barium oxide exerts the greatest influence on 
d at } agree with those of Hamilton and his coworkers [2], | increase of refractive index. The isofracts lie nearly 
and : vho have made a detailed and accurate study of parallel to the borosilicate binary , 
ander ces and densities of annealed glasses in the sys- A self-adjusting, bridge-type controller suitable for 
thy btem. Their index values are higher, in general, by a | precise temperature control of a quenching furnace 
sing } <mall but irregular amount from those of the present | was described briefly 
; 
med : apie 7. Ternary invariant points in the system BaO-B,O;-SiO 
by 
tter Ee ( 
this { 
juid ; rt j i with ome Ay . , 
uit- j BaO BLO 310 t 
s of 
: ‘ 
slo j BaO-2B203, BaO-4 B20), 810 Eut 3s 3. N ; , 
: t BaO- B2Os, BaO-2B20)5,3 Ba0-3 B:O;-2810 4 4 
ible 7 3BaO-ByOs, BaO. ByOy, BaO-SiO ‘ 19.4 t . 
BaO. ByOs, BaO- S10, 2Ba0.338i0 xf s ; vn 
ese $3 BaO0-3B9O;. 28102, BaO- 28102, 810 s f 0 x 
and $Ba0-3B:02-28i0>, Ba0-2B202,8i0 R 13.4 4 1.2 s 
bly }BaO-3ByO»-2810», BaO-28102,2Ba0-38i10 2 x4 21.4 ~ 
ork Be 8:0; Ba0 8103 28003810 | es | 1 seis 
ape ; BaO,3Ba0- ByO),2Ba0-3i0 % 
BaO-4ByO,, By, SiO Eu 2 4 uw) iM 
Phases in equilibrium with two liquids (K 
$ BaO-4BeO,, SiO) ‘ ‘ s s 
4 BaO-4B 05,8102 “ 4 
the 
ide gures 6 and 7 Because exact location of one or more boundaries 1 
>; termined by extrapolation at hed boundary curves to interse« t Postulated [21 
: ies given in weight percent Quintuple line joins compositio juids K and I 
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single-degree-of-freedom system and by 


degree-of-freedom systems 


of vibration 


mined with good accuracy 


1. Introduction 


With larger aireraft, higher landing speeds, and the 
ssitvy of flying in bad weather, the transient 
ations caused by severe gusts, landing impacts, 
ind similar shock loads are becoming increasingly 
mportant in the stress analysis of airplanes. The 
isual method of computing these transients is to 
superpose the response in a small number of the im- 
portant modes [1, 2].' Such computations are 
engthy and m some cases give results that are of 
juestionable accuracy [3] 
Houbolt [4] presents a numerical integration 
ethod in which the derivatives in the equations of 
motion are replaced by finite differences to permit a 
step-by-step calculation of the dynamic response of 
in elastic aircraft entering a gust. Houbolt’s 
method is adaptable to the problem of the dynamic 
response of an airplane to landing impact and was 
ised [5] to determine the deflections of the wing for 
an unsymmetrical two-point landing of a model air- 
From these deflections, the bending mo- 
ments on each wing at stations 17.5 inches from the 
ving root were computed. The computed results 
compared with experimental results. The 
agreement was good, indicating that the difference 
equation approach holds promise as a means of 
letermining the dynamic response of airplanes to 
shock loads. 
The purpose of the present paper is to determine 
the effect of the use of finite time increments and of 
he replacement of the continuous structure by a dis- 
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piline 


were 















e number of elastically connected point masses 
the accuracy and convergence of numerical inte- 
ition methods. The errors due to using too 
ise a time increment to adequately describe the 

detail in the foree-time history and due to 
proximating the initial conditions by various finite 
ference approximations are not considered in this 


s in brackets indicate the literature references at the end of this paper 
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rrors Introduced by Finite Space and Time Increments 
in Dynamic Response Computation 


Samuel Levy and Wilhelmina D. Kroll 


An investigation is made of the accuracy and stability of numerical integration methods 
when applied to the computation of the dynamic response of structures to impact loads 
The effect of finite time increments is studied both by 
carrying out 
the effect of finite space increments is studied by replacing a 
continuous beam by a discrete number of elastically connected point masses 

1) Of the methods investigated, only Houbolt’s is stable when the time increments are large 
compared with the natural periods of the system 
method, in this case, which result in the damping out of the responses in the higher modes 
All of the methods give good results when the time increment is less than about 
1/30 of the period in the highest frequeney mode 2 
be considered to be concentrated at relatively few mass points for computational purposes 
using a five mass idealization, the bending moment at the center of a uniform beam is deter- 


obtaining analytical solutions for a 
numerical integrations for manvy- 


It is found that 


Errors are introduced by Houbolt’s 


The distributed mass of a beam can 


report. A general study of errors in numerical inte- 
gration procedures, using techniques similar to those 
in this report, is given in references [6] and [7] 


2. Error Due to Finite Time Increment 


In order to determine the errors introduced by 
finite time increments in numerical integration meth- 
ods, a study is made of a single-degree-of-freedom 
system having simple loading conditions. The mo- 
tion of more complicated systems can be considered 
as being made up of the motion in several normal 
modes, each mode acting as a single-degree-of-free- 
dom system. The results of this study can be used, 
therefore, in judging the adequacy of the various 
numerical integration methods for complicated sys- 
tems as well as for simple systems. The results are 
determined in analytical form, using the calculus of 
finite differences, for convenience in judging the 
accuracy of the methods and the peculiar nature 
of the errors introduced by using them 


2.1. Basic Problem Considered 


The basic problem considered in this section is 
shown in figure 1. The mass, if initially disturbed, 
should vibrate without damping or amplitude build- 


up at a natural frequency w= ,hk/m (rad/sec For 
such a system, the equation of motion is 
d’r 
m -+kr=0 1) 
dt’ 


ln applying numerical integration methods to the 
solution of this equation, the following notation is 
used: 


At=time increment between successive steps in 
the numerical integration process 
n—=number of steps taken from t=0 to t= nt; 


nat 


as subscript, indicates value when f 
nat 


displacement when ¢ 
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Froure 3 ¢ app 
concerned, a study is made of the error resulting 
when this method of numerical integration is used 
and a finite number of lumped masses replace the 
distributed mass. This is done by comparing the 
computed responses with 3, 5, and 7 lumped masses 
for an undamped, uniform, free-free beam of mass 
m, length 2/ and bending stiffness, //, figure 2, using 
Houbolt’s method as modified in reference [5]. It is 


assumed that the beam is subjected to suddenly ap- | 


plied transverse loads of the type shown in figure 3. 
As more complicated time histories of loading can be 
considered to be made up of stepwise changes of the 
type shown in figure 3, the results will be significant 
for any type of transverse loading. All displace- 
ments are taken in a vertical direction. Translation 
and rotation of the beam are considered in addition 
to bending. 


3.1. Uniform Beam With Load Applied at Center 


A uniform beam is considered to be suddenly sub- 
jected to a constant normal force P at its center. 
The beam is idealized by considering the mass of 
the actual beam, figure 4 (a), to be concentrated at 
3-, 5-, and 7-mass points, figures 4 (b), 4 (ec), and 
4 (d), respectively. 

To evaluate the combined errors due to the use 
of finite numbers of mass points and finite time in- 
crements, three time increments are used in the 
numerical analysis for each mass distribution. 
These time increments are 


: l ; l ; 
M==,mG/12ET7; avml 1L2EJ, and -ym/12F].- 
. 7) 


m, EF 


s 





(a) Actual beam 


™m/4 


O 


(b) Three mass system 


m/s m/4 m/4 m/4 m/8 
Ow-Orw 
(c) Five mass system 
m/i2 m/6 m/6 m/6 m/6 m/6 m/i2 
(d) Seven mass system 
Fioure 4 Distributions of masse considered for beam subie 


to m pact at center 


The computations are carried out as described in 
reference [5} and as shown for a typical case in the 
appendix. In this procedure, the initial conditions 
are taken into account by considering the displace- 
ments zero for three time increments prior to the 
initial instant of time. The flexibility is taken inte 
account by the use of influence coefficients. Matrix 
algebra is used to simultaneously compute the dis- 
placements at all the mass points at a particular in- 
crement of time from the given forces acting at that 
time and from the known displacements at the pre- 
ceding three increments in time. The bending 
moment at the center is then computed from th 
inverse of the influence coefficient matrix, the dis- 
placements at the various mass points, and the lever 
arms from these mass points to the center of th 
beam. The bending moment, ./, at the center of 
the beam, positive when concave away from the side 
of the beam being struck, is plotted on a dimension- 
less basis in figures 5, 6, and 7 for the 3-, 5-, and 7- 
mass approximations, respectively. It is noted that, 
as the time interval used in the step-by-step nu- 
merical integration is shortened, the maximum abso- 
lute value of the bending moment is increased. As 
is to be expected from the analysis of a single de- 
gree-of-freedom system, the errors of the numerical 
method cause a damping out of the oscillatory re- 
sponse and, as the time increment is decreased, a 
decrease in the period of the oscillatory response. 

For each of figures 5 to 7, the mass distribution is 
kept constant, and the increment of time used in 
the numerical integration is varied to show the effect 
of time increment. In each of figures 8 to 10, the 
time increment is kept constant, and the responses 
for the different mass distributions are plotted to 
bring out more clearly the effect of increasing th: 
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number of masses. The maximum absolute value 
of the bending moment at the center differs by not 
more than 4 percent for the 5- and 7-mass systems, 
regardless of time increment. There is some differ- 
ence in the period of the response for the shorter 
time increments for the 5- and 7-mass systems 


3.2. Beam With Impact at One Tip 


Because it was believed that an impact at a tip 
mass might result in a more severe loading condition, 
an investigation is made for the case of a load P 
applied at the tip of the beam shown in figure 4 (a 
using the idealizations shown in figures 4 (b) and 
4 (c). The bending moments at the center for the 
3- and 5-mass idealizations are shown in figures 11 
and 12, respectively. It is noted again, as in case 
3.1, that as the time interval used in the step-by- 
step numerical integration is shortened, the maximum 
absolute value of the bending moment at the center 
is increased. As is to be expected from the analysis 
of a single-degree-of-freedom system, the errors of 
the numerical method cause a damping out of the 
oscillations, which is particularly evident for the 
longest time increment, and a decrease in the period 
of the oscillatory response as the time increment 
is decreased. 

In figures 13 to 15, the time increment is kept 
constant, and the responses for idealization of the 


beam by 3 and 5 masses are plotted to bring out 
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fect of increasing the number of masses. It 
be seen that the period of the initial vibratory 
ise decreases with an increase in the number 
masses and that the damping increases. The 
aximum bending moment at the center is relatively 


lait cted. 


13.3. Pair of Beams Representing a Fuselage and 
Wing With Impact at Both Wing Tips 


\n airplane can be approximated by two uniform 
eams at right angles to each other, one beam repre- 
senting the wing, the other the fuselage. Such a 
Each beam 


i system, figure 16 (a), is investigated. 
s idealized as 5 masses, as sbown in figure 16 (b). 


suddenly applied constant loads, ?, are consid- 
bered to act at each wing tip. The bending moments 
} at the center of both the ‘‘fuselage”’ and ‘‘wing”’ beams 
ire computed by using the same three time incre- 
ents in the numerical integration as were previ- 

used. The bending moment at the center of 

the struck beam (wing) is shown in figure 17 and 
that for the other beam (fuselage) is shown in figure 
8. The agreement between the bending moment 
ratios for the three values of time increment is only 
far even for the two shortest time increments. As 
might have been expected from the analysis given 
of the single-degree-of-freedom system, the error 
due to using the longest time increment is evident 
the marked damping of the oscillatory response. 


3.4. Discussion 


lt is estimated that a time increment of about 1/30 
{ the period in the fundamental mode is necessary 
lor good accuracy in numerical-integration methods. 
his value is based on several factors. 

The solution given in eq (26) corresponds to 
about 20 time increments per fundamental period. 
In this ease the frequency coefficient, 0.961, is in 
error by about 4 percent. 

The convergence indicated in the figures is 

lite complete for the shortest time increment, | 
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which corresponds to 1/13.74 of the fundamental 
symmetrical mode period and only 1/5.00 of the 
lowest antisvmmetrical mode period. A shorter 
time increment than this is certainly needed to give 
the response of the higher modes 

(3) The recommended time increment is likely to 
be used for computing responses with slower con- 
vergence than the bending moment at the center of 
the beam and for structures with slower convergence 
than that of a uniform beam 

(4) It is likely that computations by numerical 
methods will be carried out by automatic computing 
machines where the advantage of additional accuracy 
from a smaller time increment can be had with less 
penalty than when hand computing is used 

It is noted from figures 8, 9, and 10 that, for each 
time increment considered, the response to an impact 
at the center is essentially the same for the 5- and 
7-mass idealizations of the beam but considerably 
higher for the 3-mass system. For impacts at one 
end of the beam, figures 13, 14, and 15, the response 
shows less change in going from the 3- to the 5-mass 
idealizations than is the case for impact at the center 
It would appear that a uniform beam should be 
approximated by a 5-mass idealization if errors in 
the bending moment at the center are to be small 

In table 1, the maximum absolute values of 
bending moment ratios are given. For the beam 
with the impact at the center, the absolute value of 
the bending moment increased as the time inter- 
val for integration was decreased, but decreased 
as the number of masses used for the idealization of 
the distributed mass of the beam was increased 
From this, it may well be that the errors caused by 
using a reasonable value of time interval for the 
numerical integration will be offset to some extent 
by the errors caused by using a finite number of 
masses. For the impact at a tip mass, the increase 
in bending moment with decrease in time interval is 
appreciable, but the decrease in bending moment with 
increase in number of masses is much less than for the 
case of impact at the center 
wal we moment 


absol ute of bhendi ng 


with force P 


TABLE 1 Varimum 
center of uniform beam struck 


Impact at center 


Impact 


Vy mi i2 
1/3y mi 


Va\m 








It is evident from figure 5 that the period of the 
response approaches the value for the exact solution 
for a beam idealized by 3 masses as the time interval 
in the integration is shortened. This same decrease 
in period with time increment is noted in figures 6 
and 7 for the 5- and 7- mass idealization of the beam. 


The period of the response also decreases as the num- | 


ber of masses is increased, figures 7, 8, and 9 
latter result, at least qualitatively, verifies that ob- 
tained by Dunean [9]. He has obtained an exact 
solution showing that the periods of the computed 
modes of vibration of a uniform beam decrease as 
the number of masses chosen to represent the beam 
is increased, and that the magnitudes of the errors 
in the computed periods vary inversely as the square 
of the number of masses. 


4. Conclusions 


This | 


From the analyses of the responses of single-degree- | 


of-freedom and many-degree-of-freedom systems to 
a suddenly applied force, it may be concluded that: 
a) Numerical integration methods give results 


that closely approximate the exact solution only if | 


the increment of time between successive steps in 
the integration is small compared to all the natural 
periods of vibration of the system 
to those investigated where bending moment near 
the center is desired, a time increment of about Yo 
of the fundamental period of the system will give 
good results 

b) Of the methods investigated, the Houbolt 
method is the only one that gives convergent results 
for large time increments. In this method, errors 
introduced by large time increments result in the 
damping out of oscillatory response 

c) The distributed mass of the beam can be con- 
sidered to be concentrated at relatively few mass 
points for computational purposes. In cases similar 
to those investigated, approximating a uniform beam 
bv a 5-mass idealization gives good results. 


5. Appendix 


The detailed method of determining the response 
of a uniform beam to an impact load by numerical 
integration is given in this appendix for the case of 
a beam subjected to an impact at its tip. The case 
considered is that where the time increment is 
At=ymP/12ET. The responses for all the other 
cases presented in this report are obtained in a 
similar manner 
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" ° | 
In cases similar | 


| 


The mass of the uniform beam is considered | 
distributed at five stations along its length, figu 
The force, P, is applied at the tip, station 5 
used in this report, an influence coefficient, 6, , is 
displacement at station r due to a unit load at st; 
s with the center of the beam, station 3, fixed. 
influence coefficients are given in table 2. Deng 
the displacement positive upward at station r b 
the force at station r by F,, and the rotation at sta 
3 by a, the displacements of the beam at statior 
2,4, and 5are 


d 


d,- 


d;+ 





d, d, T 


Substituting the values of the influence coefficient: 
from table 2 into eq (Al), and solving for the forces 


, - ) 
sence coefficients 6, , for 5-mass idealized be 


shown in figure 19 


TABLE 2 Infl 


in terms of the displacements at each station and \ 
the rotation at station 3. 


isk] 


A= "ap 


( 2d, 5d, +3d, ~ )> 


48k] 
7f 


(—5d,+ 16d,—11d;+3ad), 
48ET 


7 Sd’ 


al 
) | 


lid, + 16d, 3a), 


3d;—5d,+-2d;4 





48k] 
cp ( 


From the condition that the sum of the forces equa 
zero, 


and substituting from eq (A2), the value of F; is 


Fy AO (3dy—11ds + 16d,— 11d,4 3d (Ad 





the requirement that the net moment equal 
» Ae ; , f3l 
0)+F:(5)+FiO+F. (5 )+Fs(20=0, (A5 


(A4), 


sing eq (A2) and 


18k] 


5d. 
7/ ae, 


1.875d 4.2 


——- $.25d,4+-11.50d, 


48k] 


7] 11.00d.- 


(3.00d, 


418k 
7/ 
48K] 


F =) (0.125d, 


F, 0.75d,+-4.50d, 
0.75d,4 


(he forces are the sum of the inertia loads due to 
he acceleration of the mass of the beam and the 


> applied forces, 


{| P 


4 





acceleration of 


and P, applied 


here m is mass of uniform beam; d,, 
wam at station r(r=1,2, ... 5); 
onstant force at station 5. 
i 
i 


Using the difference eq , acceleration at station 


can be expressed as 


(AQ) 


quating the first of eq (A7) and (AS) and substi- 
g for (d,), from eq (A9), we obtain 


[2(d,), 5(di)s ait 4(d;), 24t~ di)s sacl 


inf 
i (1.875d, —4.25d,.+-3.00d,—0.75d,+0.125d;5];. 


(A10) 


3.00d 


16.00d 


11.00d,4+-11.50d, 


3.00d 


A2), 


Substituting the value of a from eq (A6) into eq 
and rewriting eq (A4 


0.75d, T 0 125d 


11.00d, ; 4.50d, 0.75d-). 


11.00d,4-3.00d,), 
4.25d;), 


$.25d,+-1.875d;) 





Similar relations are obtained by using the remain- 
ing eq (A7) and (AS). 

Multiplying eq (A10) and the corresponding equa- 
tions for the other stations through by 3/7/P/ 
express the deflections as a dimensionless ratio, and 
transposing terms, the equations can be written in a 
form which, for station 1, reads 


7Tml 
192A it 


3d,E] 
PR / 


to 


3d,E1 


3d,E]1 
)- Pr 


(1.8754 Pp 


4.25 ( 


ad 3d,k]1 
0 75 ( Pl ) 


)4 
3.00 ( 


3d,E1 
pp) 


) 125 ( 


where 


3d, EI 
PP), 7 


) | Al2) 
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Equation (A111) and the equations similar to it can be 
expressed in matrix form as 


7Tml a ( 
38442 7ae |” 


3d,E1 


{ 3d,E] Pp 


Pl 


\ 3d,E] ) 


Al) pp ¢, 


(A13) 


where the matrix [A] is given in table 3. Solving eq 


(A13) for the deflections gives 


3d,E1) 


( 7 
) pp (oll 








TABLE 3 ase in equation (A138 


where [A]~' is obtained as the inverse of [A] with | where 


At=,mP 12k]. The matrix [A]~' is given in table 4 
The first of this set of equations in conventional form 


Is 
3E] 


~ 3d,E1 ) 40.250 3d,k 1 ) \ 


Pl Pl 
3d,E1 . 3d,E1 
PI ) +0.250 ( ph ) + 
3d,E] 

py), 


2m D,—0 125 ( 
3d,E1 
(“pp ) 


0.113294(7 0.113732(.J 0.237485(J 


1.344937(J,), +0.589567(.J7, 


0) 250 ( 


where, from eq (A12), Q 125 ( 


3d Kk] ) 0. X7500 sd kl ) 

PI : ving Pi oo The other equation is obtained by substituting 
EI (3d, ED)/PP),, (Bd, /PP),, . 2. [(8d,ED)/P 
a ) Al6) | in eq (A19), their equivalents in terms of the J's 
Pl Al4 This gave 


| 09375 ( 
+ 0.21875 ( 


Similar equations can be obtamed for J at the other 3k] 
ae 28: 375( J), +0.28571362(J; 


stations PI 
t ().28571362(/, 


\2 


The small differences in the coefficients in eq (A2 
result from rounding-off error. For a good e| 
indication, the coefficients, as given, are need 
It is assumed that, prior to the application of the The values of D from = A19) and \20) are gi 
force, ?, the displacements at all stations are zero in columns (7) and (8) of table 5 
The computations are tabulated as shown in tables The J’s are checked by comparing their sums 
. eo . ' letermined from two different formulas. The f{ 
5 and 6. The J’s are computed for time f, table 6, | ' 
using eq (A16) for (J,), and similar equations for 
the other J’s. The d’s are computed for time tf, 
table 5, using eq (A15) for (3d,47/PP), and similar (>> J,) af P wer L(T) \ 
equations for the other d’s. From the d’s at time f¢, oa 7 ve 4 
the J’s at time t+ At are computed and, from these, 
the d’s at time t+At. Thus a time history of the | The second formula is obtained by substituting in 
displacement ratios of the uniform beam due to the | (A21) the value of (J;), in eq (A16) and correspo 
suddenly applied constant force is obtained. ing values for the other J’s. Then, making use of 
Since the step-by-step numerical procedure is | (A419). 
dependent on previously computed values, it is 
desirable to have a check column which will indicate 
errors made in the computation. The deflections at ( 
each station were checked by computing the dis- 
placement, J), of the center of gravity of the beam 1.75(D) 10.4375 \2 
from two different formulas. If these two values a ted Seaman a 
of )) agreed, the computation of the deflections was | The values computed from eq (A21) and (A22) a 
assumed correct. The first equation for D, is: given in columns (7) and (8), respectively, of tal 
3K] , 3d.E] 6. The good agreement between the check columm 
Pi dD, = ( m,: PR )» (A17) | in tables 5 and 6 indicate that the computation * 
r=l free of numerical error. 
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The forces are determined from the displacements 
of the masses at the five stations on the beam by use 
of eq (A7). The forces, multiplied by thier distances 
from the center of the beam determine the bending 
moment at the center of the beam 
bending moment at the center by 17, 


l 


») 


M=F/+ F; (A23 


2 50d,4 
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Values of the bending moment ratio computed 
from eq (A25) are given in column 9, table 5 


They 
are plotted in figure 12 for At= ymP/12Fk7 
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puting the responses of the uniform beams to impact 
loads 
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